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ABSTRACT 

Davidson,  John  Leitch.  Ph.D.,  Purdue  University, 

August  19  73.   The  Effect  of  Quasi-Preconsolidation  on 

Compression  of_  Clay  Soils .  Major  Professor: 
Dr.  William  II.  PerlofT^ 

The  research  was  initiated  to  predict  pore  pressures 
and  displacements  during  primary  consolidation  of  a  clay 
soil  which  exhibits  a  quasi-preconsolidation  effect.   A  soil 
model  representing  this  behavior  was  formulated  from  labora- 
tory testing,  and  was  incorporated  into  a  finite  element 
program  to  yield  solutions  to  realistic  boundary  values 
problems . 

In  the  laboratory  study,  drained  triaxial  tests  were 
performed  on  soil  specimens  in  which  a  quasi-preconsolida- 
tion effect  had  been  induced  by  maintaining  them  under  a 
constant  hydrostatic  stress  for  a  period  of  time.   Resulting 
stress-strain  plots  showed  an  essentially  bi-linear  behavior, 
with  an  initially  high  and  subsequently  lower  stress-strain 
modulus.   The  break  point  represents  yielding  of  the  bonds 
formed  during  the  period  of  constant  loading,  and  is  equiva- 
lent to  the  phenomenon  observed  by  Leonards  and  Ramiah  during 
consolidation  testing  and  designated  by  them  as  the  quasi- 
preconsolidation  pressure  point.   From  different  stress 
path  tests  on  specimens  consolidated  under  different 
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hydrostatic  pressures,  a  three-dimensional  yield  surface  was 
formulated  in  terms  of  the  hydrostatic  stress  and  the  test 
effective  octahedral  stress  increments. 

The  finite  element  solution  to  the  consolidation  prob- 
lem, based  on  the  Biot  theory  for  an  elastic  soil,  calculated 
stresses,  pore  pressures  and  displacements  within  the  soil 
after  each  increment  of  pore  pressure  dissipation.   A  bi- 
linear model  with  associated  yield  surface,  derived  from 
laboratory  test  results,  was  introduced  into  the  finite 
element  program  and  solutions  were  obtained  for  one-dimen- 
sional and  two-dimensional  plane  strain  consolidation  prob- 
lems. 

Solutions  were  obtained  using  the  Purdue  University 
Control  Data  Corporation  6500  Computer  System  and  employing 
a  Fortran  IV  source  deck. 

Results  from  the  finite  element  programs  incorporating 
the  yield  model  indicate  the  following  soil  behavior. 

1.  Prior  to  yielding  settlements  are  small,  pore 
pressures  are  small  and  dissipation  rapid. 

2.  Upon  yielding  pore  pressures  increase  due  to  the 
transfer  of  load  from  the  soil  skeleton  to  the  pore 
water.   Settlements  are  large  and  pore  pressures 
dissipate  slowly. 

3.  Such  behavior  is  in  agreement  with  recently  pub- 
lished field  instrumentation  studies. 


1.  INTRODUCTION 

In  order  to  satisfactorily  predict  stresses  and  deforma- 
tions in  a  soil  mass  which  will  act  as  a  foundation  for  con- 
struction, a  number  of  interrelated  phases  must  be  considered. 
A  subsurface  investigation  of  the  site  yields  both  samples 
and  a  delineation  of  the  insitu  soil  conditions.   Laboratory 
and  field  tests  are  then  performed  on  representative  samples 
in  order  that  soil  parameters  may  be  defined,  and  these  and 
their  variations  are  used  in  the  formation  of  a  descriptive 
model  or  material  constitutive  law.   A  solution  procedure 
incorporating  the  model  is  employed  to  calculate  deformations 
and  stresses  for  the  specific  boundary  value  problem.   Finally, 
measurements  from  field  instrumentation  during  construction 
are  compared  with  predicted  values  and,  if  necessary,  lead  to 
alteration  of  design  or  construction  procedures. 

Much  research  in  recent  years  has  been  directed  towards 
better  formulation  of  the  descriptive  model  and  towards 
development  of  better  and  more  versatile  solution  procedures 
for  such  models.   Traditionally,  earth  masses  have  been 
modelled  as  homogeneous  isotropic  materials,  rigid  and  per- 
fectly plastic  for  stability  analyses,  and  linear  elastic 
for  stress  distribution  and  settlement  analyses. 

In  current  research,  soils  are  being  more  realistically 
modelled  as  inhomogeneous  anisotropic  materials,  characterized 


by  non-linear  elasticity,  viscoelasticity ,  strain  hardening 
and  softening  plasticity  and  combinations  thereof. 

A  major  problem  exists,  however,  in  relating  models 
derived  in  the  laboratory  on  samples  deformed  under  simple 
stress  conditions  to  more  complex  field  problems.   Exact 
solution  procedures  are  often  capable  of  considering  only 
highly  idealized  soil  constitutive  behavior  and  simple 
boundary  conditions. 

In  the  past  decade,  advances  in  the  field  of  computer 
sciences  and  the  ready  availability  of  large  digital  compu- 
ters have  transformed  much  of  the  emphasis  on  solution  pro- 
cedures from  exact  to  numerical  methods.   With  the  aid  of  the 
computer,  procedures  such  as  the  finite  difference  and  the 
more  recent  and  versatile  finite  element  method  become  effici- 
ent and  accurate  methods  of  solving  field  problems  with  com- 
plicated boundary  conditions  and  more  realistic  soil  behavior. 

However,  models  are  still  by  necessity  simplifications  of 
reality,  and  must  be  chosen  to  reflect  those  characteristics 
deemed  most  significant  for  the  problem  at  hand.   This  research 
represents  a  study  of  one  such  behavior,  namely  the  quasi- 
preconsolidation  effect.   Based  on  laboratory  results,  a 
stress-deformation  model  considering  this  behavior  is 
formulated.   Incorporation  of  the  model  into  a  finite  ele- 
ment program  allows  realistic  boundary  value  problems  to  be 
solved. 


2.   REVIEW  OF  LITERATURE 

2.1   Introduction 

Karl  Terzaghi's  formulation  of  the  one-dimensional 
consolidation  theory  is  generally  considered  the  beginning 
of  modern  soil  mechanics.   And  yet,  the  problem  of  predicting 
the  magnitude  and  rate  of  deformation  of  soft  fine  grained 
soils  has  remained  one  of  the  main  concerns  and  areas  of 
research,  both  theoretical  and  practical,  in  geotechnical 
fields.   This  is  especially  so  now,  when  "unsatisfactory" 
subsurface  conditions  often  play  a  minor  role  in  determining 
the  location  and  size  of  engineering  structures  such  as  large 
buildings  and  highway  embankments.   Foundation  engineers, 
more  so  now  than  ever,  are  required  to  use  the  material 
available  at  the  site  available. 

The  process  of  consolidation  can  be  viewed  as  a  transfer 
of  load  from  pore  water  of  low  compressibility  to  a  more  com- 
pressible soil  structure.   Loads  applied  at  the  soil  bound- 
aries create  pressures  in  the  pore  fluid  in  excess  of  hydro- 
static, and  gradients  in  this  excess  pressure  induce  flow 
towards  drainage  zones.   This  physical  process  whereby  pore 
water  is  slowly  squeezed  out  of  the  clay  stratum  due  to 
external  loads  has  long  been  recognized.   In  1809  Telford  in 
Scotland  preloaded  a  thick  soft  clay  layer  for  nine  months 


before  constructing  a  canal  lock.   In  the  United  States, 
papers  by  W.  S.  Smith  in  1892  and  E.  C.  Shankland  in  1896 
described  building  settlements  in  Chicago  due  to  consolida- 
tion.  However  a  theoretical  explanation  and  reliable 
methods  for  estimating  consolidation  settlement  and  time 
rate  did  not  exist  until  Terzaghi  first  introduced  the 
principle  of  effective  stress  in  1921  and  his  classical 
theory  of  consolidation  two  years  later  (Terzaghi,  1923). 

The  one-dimensional  equation  derived  by  Terzaghi, 
governing  the  dissipation  of  excess  pore  pressure  with  depth 
and  time  is 

a  linear  partial  differential  equation  of  parabolic  form 
where 

u(z,t)  -  excess  (above  hydrostatic)  pore  water  pressure 
t  -  time 

z  -  depth 

v 

c  B  r; -  coefficient  of  consolidation 

v   Y  m 
w  v 

k  -  permeability 
Y  -  unit  weight  of  water 
m  -  coefficient  of  volume  compressibility 

This  equation  is  of  the  same  form  as  that  for  the  non- 
stationary  one-dimensional  flow  of  heat  through  isotropic 
bodies  and  consequently  a  complete  methodology  and  a  large 
number  of  solutions  became  immediately  available.   Carslaw 


(1921)  in  his  book  "Introduction  to  the  Mathematical  Theory 
of  the  Conduction  of  Heat  in  Solids"  outlines  solutions  to 
the  heat  flow  equation  for  many  different  boundary  and  ini- 
tial conditions.   A  closed  form  solution  based  on  Fourier 
mathematics  is  given  by  Terzaghi  (1941)  for  boundary  and 
initial  conditions  corresponding  to  the  laboratory  consoli- 
dation test. 

In  order  to  reduce  the  complex  consolidation  problem  to 
the  above  simple  equation  it  was  necessary  for  Terzaghi  to 
make  a  number  of  linearizing  assumptions.   These  greatly 
limit  the  generality  of  the  problems  which  can  be  solved  and 
researchers  have  modified  them  in  deriving  solutions  to  prob- 
lems with  particular  soils  and  boundary  conditions  which 
depart  significantly  from  those  assumed  by  Terzaghi. 

There  are  three  basic  assumptions  inherent  in  the  deri- 
vation of  the  Terzaghi  consolidation  equation.   First  is  the 
premise  that  the  consolidation  settlements  are  related  to 
and  dependent  solely  on  the  dissipation  of  the  excess  pore 
water  pressure  developed  by  the  loading;  that  is,  a  hydro- 
dynamic  process.   Secondly,  it  is  assumed  that  the  soil 
deposit  is  both  isotropic  and  homogeneous.   No  spatial  vari- 
ation of  material  parameters  is  possible,  and  both  deforma- 
tion and  drainage  are  assumed  to  be  one-dimensional.   Final- 
ly, in  order  to  achieve  the  simple  resulting  differential 
equation  it  was  necessary  for  Terzaghi  to  assume  a  constant 
coefficient  of  permeability  and  a  linear  stress-strain,  or 
void  ratio-pressure  relationship. 


The  first  assumption  regarding  the  hydrodynamic  effect 
has  been  generally  accepted  by  subsequent  researchers  and 
efforts  have  been  concentrated  on  relaxing  the  homogeneity 
and  linearity  assumptions.   A  brief  review  of  such  research 
is  given  in  the  following  two  sections. 

2. 2  Spatial  Variation  of  Soil  Parameters 
A  number  of  solution  procedures  has  been  introduced  to 
take  into  account  the  non  homogeneous  nature  of  virtually 
•all  soil  deposits.   Parameter  variation  with  depth  can  often 
be  approximated  by  a  number  of  essentially  homogeneous  hori- 
zontal strata:   the  n-layer  problem.   This  has  been  solved 
in  closed  form  (Gray,  1945),  by  averaging  coefficients 
(Terzaghi,  1940),  and  by  a  transformation  method  (Glick, 
1945) .   The  finite  difference  solution  procedure  (Forsythe 
and  Wasow,  1960)  has  been  adapted  for  the  one-dimensional 
consolidation  problem  and  used  by  several  researchers 
including  Abbott  (I960)  ,  Raymond  (1966)  ,  and  Davis  and  Lee 
(1969). 

Other  authors  have  published  solutions  to  problems  in 
which  particular  soil  parameters  were  permitted  to  vary  in 
some  manner  with  depth.   Edelmann  (1953)  obtained  an  analyti- 
cal solution  to  the  problem  in  which  the  modulus  of  elastic- 
ity increased  with  depth.   Schiffman  and  Gibson  (1964)  con- 
sidered the  case  in  which  the  permeability  and  the  coeffici- 
ent of  compressibility  varied  with  depth,  while  Martins 
(1965)  and  Schiffman  (1958)  permitted  only  the  permeability 
to  vary. 


Two-  and  three-dimensional  consolidation  theories  have 
been  formulated,  based  on  the  work  of  Biot  (1941),  and  solu- 
tions derived.   These  have  however,  by  necessity,  been  solu- 
tions to  problems  with  simple  boundary  configurations,  simple 
loading  conditions  and  ideal  constitutive  laws.   No  solution 
procedure  existed  which  could  take  into  account  such  complex 
conditions  until  the  Finite  Element  Method  was  introduced, 
(Turner  e_t  al_.  ,  1956;  Clough,  1960;  Wilson,  1965).   Using 
this  versatile  numerical  tool,  many  solutions  to  specific 
soil  engineering  problems  have  since  been  published.   A 
review  of  many  of  these  papers  is  given  by  Radhakrishnan  and 
Reese  (1970)  ,  and  illustrates  the  versatility  of  the  Finite 
Element  Method  in  areas  of  non-linear  constitutive  behavior, 
random  spatial  variation  of  soil  parameters,  and  complex 
boundary  and  loading  conditions. 

2 . 3  The  Linearity  Assumption 
Formulation  of  the  consolidation  equation  by  consider- 
ing flow  into  and  out  of  a  soil  domain  results  in  the  follow- 
ing equation 
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where 

e  -  void  ratio 


e   -  initial  void  ratio 
o 


and  other  terms  are  as  defined  for  equation  (2.1). 


In  order  to  solve  for  the  excess  pore  pressure,  u(z,t),  it 
is  necessary  to  assume  relationships  between  permeability 
and  void  ratio  or  permeability  and  effective  stress,  and 
between  void  ratio  and  excess  pore  pressure  or  void  ratio 
and  effective  stress. 

Terzaghi  in  his  fundamental  derivation  assumed  a 
constant  coefficient  of  permeability  and  the  linear  void 
ratio-pressure  relationship, 
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Schiffman  (1965)  assumed  a  linear  void  ratio-pressure 
relationship  but  permitted  the  permeability  to  decrease 
either  step  by  step  with  time,  or  as  an  exponential  function 
of  the  excess  pore  pressure. 

Davis  and  Raymond  (1965)  assumed  a  linear  void  ratio- 
logarithm  of  pressure  relationship  and  varied  the  permeabil- 
ity such  that  the  coefficient  of  consolidation  c  remained 
constant. 

Barden  and  Berry  (1965)  ,  Raymond  (1969)  ,  and  Poskitt 
(1969)  all  assumed  linear  void  ratio-logarithm  of  pressure 
and  linear  void  ratio-logarithm  of  permeability  relationships 

Gibson  e_t  al.  ,  (1967)  formulated  a  consolidation  equa- 
tion in  which  both  the  permeability  and  the  effective  stress 
were  general  functions  of  the  void  ratio.   This  reduced  to  a 
form  similar  to  Terzaghi 's  equation  but  with  a  more  complex 
coefficient  cp.   Solutions  were  obtained  assuming  cp  constant 
and  linearly  related  to  the  void  ratio. 


Results  from  the  above  theories  have  been  compared  with 
the  simple  linear  Terzaghi  formulation,  and  in  certain  cases 
do  show  some  variation.   However,  it  was  noted  by  Barden  and 
Berry  (1965)  that  at  the  field  scale  the  simple  linear 
Terzaghi  theory  may  generally  be  applied  to  the  majority  of 
clays  exhibiting  no  unusual  properties  such  as  partial 
saturation  or  high  compressibility. 

Raymond  (1966)  performed  consolidation  tests  on  two 
different  soils  in  an  apparatus  consisting  of  either  two  or 
four  oedometer  cells  with  their  drainage  lines  connected  in 
series.   He  observed  an  apparent  overconsolidation  behavior 
during  consolidation  after  a  load  was  applied.   He  formulated 
a  one-dimensional  finite  difference  program  capable  of 
describing  the  consolidation  process  during  a  pressure  incre- 
ment in  which  the  preconsolidation  pressure  was  exceeded. 
Assumptions  were  made  of  a  bilinear  void  ratio-pressure 
relationship  with  the  break  at  the  preconsolidation  pressure, 
and  of  constant  but  different  coefficients  of  permeability 
above  and  below  this  pressure.   The  consolidation  process 
then  consisted  of  an  internal  boundary  moving  through  the 
initially  overconsolidated  soil.  As  this  boundary  passed 
the  soil  behavior  changed  from  overconsolidated  to  normally 
consolidated.   Theoretical  results  showed  that  consideration 
of  even»-s lightly  overconsolidated  behavior  resulted  in  sig- 
nificant differences  from  the  Terzaghi  solution,  indicating 
a  rapid  dissipation  of  pore-water  pressure  until  a  critical 
point  was  reached  after  which  dissipation  was  much  slower. 
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Many  clays  which  from  their  geologic  history  would  be 
classified  as  normally  consolidated  do  not  behave  in  the 
expected  simple  linear  void  ratio-logarithmic  pressure 
manner  when  initially  loaded,  but  rather  show  an  apparent 
overconsolidation  behavior  of  the  type  observed  by  Raymond. 
In  fact  it  appears  that  most  clay  soils  if  permitted  to 
remain  under  a  constant  stress  state  for  even  a  short  period 
of  time  will,  in  some  manner,  develop  interparticle  bonding 
which  can  support  some  additional  loading  with  very  little 
deformation.   If  however,  the  applied  pressure  exceeds  some 
critical  value  the  comparatively  rigid  bonds  are  susceptible 
to  sudden  breakdown.   This  behavior  results  in  a  sharp  slope 
increase  on  the  e-log  p  curve  at  a  pressure  which,  by  analogy 
with  the  effect  of  precompression,  has  been  termed  the 
"quasi-preconsolidation  pressure,  pca"  by  Leonards  and  his 
co-workers  at  Purdue  University.   This  phenomenon  is  dis- 
cussed in  detail  in  the  following  section. 

2.4  The  Quasi-preconsolidation  Effect 
A  normally  consolidated  soil  deposit  is  defined  as  one 
in  which  the  existing  effective  overburden  pressure  p0  is 
equal  to  the  maximum  effective  pressure  that  the  soil  has 
ever  experienced  pc.   Figure  2.1  shows  a  schematic  void 
ratio-vertical  pressure  plot  of  normally  consolidated  clay 
behavior.   Line  AB  represents  the  geologic  sedimentation  of 
the  soil  and  BC  the  projected  behavior,  should  further 
loading  be  applied.   Point  B  denotes  the  present  insitu 
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FIGURE     2-1         SCHEMATIC     VOID     RATIO  -  LOGARITHM 
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state,  with  void  ratio  e   and  vertical  pressure  equal  to  the 
overburden  pressure  p  .   Curve  DC  represents  the  results  of 
a  consolidation  test  on  an  "undisturbed"  sample  of  the  clay 
and  consists  of  a  recompression  curve  which  bends  into  the 
virgin  compression  curve  at  a  pressure  p  equal  to  the  over- 
burden pressure. 

Since  BC  is  approximately  linear  on  the  semi-logarithmic 
plot  settlement  calculations  are  made  using  the  equation 

o   C       p  +  Ap 

«  -  I  TTJr  i°*  ^  »» 


Pq 


where 

p  is  the  effective  overburden  pressure 

Ap  is  the  vertical  load  increment 

Az  is  the  thickness  of  an  individual  clay  layer 

Cc   is  the  compression  index  (slope  of  the  virgin  com- 
pression curve) 

Bjerrum  (1964)  published  a  report  in  which  he  compared 
observed  settlements  of  structures  and  embankments  on  soft 
clay  with  computations  made  using  the  above  equation.   The 
results  showed  that  for  heavy  structures  with  large  settle- 
ments good  agreement  was  obtained,  but  for  light  structures 
the  calculations  greatly  overpredicted  the  actual  field 
settlements.   The  explanation  of  this  discrepancy  lies  in 
the  failure  of  the  analysis  to  take  into  account  the  quasi- 
preconsolidation  effect  which  exists  in  most  clay  soils. 
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This  effect  manifests  itself  on  the  void  ratio-log  pressure 
plot  (Figure  2.2)  as  an  almost  horizontal,  low  compressibil- 
ity line  from  p  to  p   when  the  effect  is  destroyed.   The 
curve  then  reverts  to  the  approximately  linear  virgin  com- 
pression line.   It  is  this  latter  fact  which  accounts  for 
the  large  structure  settlement  agreement  reported  by  Bjerrum, 
since  loading  was  of  sufficient  magnitude  to  destroy  and 

exceed  any  p   effect.  However,  for  a  small  load  increment 
cq 

failure  to  consider  the  initial  stiffness  results  in  large 
errors.   These  two  loading  cases  are  shown  schematically  in 
Figure  2.3. 

Extensive  research  on  the  quasi-preconsolidation 
phenomenon  has  been  carried  out  by  Leonards  and  others  at 
Purdue  University.   A  program  of  one-dimensional  consolida- 
tion tests  on  two  remolded  soils,  a  limestone  residual  clay 
and  a  silty  clay  of  glacial  origin,  was  performed  and 
reported  by  Leonards  and  Ramiah  (1960).   During  testing,  a 
single  pressure  increment  was  maintained  for  periods  of  time 
up  to  100  days.   In  some  tests  secondary  compression  during 
this  period  was  permitted,  while  in  other  tests  after  a  24 
hour  period  of  primary  consolidation  further  volume  change 
was  prevented.   In  all  cases,  whether  secondary  compression 
was  permitted  or  not,  continued  loading  after  the  rest 
period,  using  small  pressure  increments,  produced  a  quasi- 
preconsolidation  effect  considerably  larger  than  the  maximum 
previous  consolidation  pressure.   Figure  2.4  from  Ramiah 
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Overburden    Pressure 

Quasi-Preconsolidation    Pressure 


Load    remains  constant  for 
some    long    period   of    time 


Virgin   Compression- 
Curve 


Logarithm    (Pressure) 


FIGURE    2-2     THE     QUASI-PRECONSOLIDATION 
EFFECT 
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(1959)  illustrates  this  phenomenon  which  was  observed  in  all 
tests  on  both  clays  and  at  varying  moulding  water  contents. 
Further  tests  were  carried  out  in  which  samples  were  pre- 
pared as  before  except  that  after  the  rest  period  the  speci- 
mens were  unloaded  and  "sampled"  by  rapid  penetration  of  a 
floating  ring  consolidometer  with  sharpened  edge.   Consolida- 
tion tests  were  then  performed  and  the  quasi-preconsolidation 
effect  again  observed. 

Further  experimentation  was  performed  on  artificially 
sedimented  samples  of  the  limestone  residual  clay  (Leonards 
and  Altschaeffl,  1964).   Samples  were  prepared  and  loaded  in 
lucite  tubes  at  rates  approaching  those  encountered  in  the 
formation  of  river  deltas.   Different  combinations  and  rates 
of  loading,  unloading,  reloading  and  rest  periods  were 
easily  accomplished  in  the  especially  designed  sedimentation 
apparatus  (Altschaeffl,  1960).   Figure  2.5  shows  the  void 
ratio-log  pressure  plot  for  one  test  which  was  loaded  to 
0.512  kg/cm2  and  then  allowed  to  remain  under  this  constant 
pressure  for  100  days.  The  specimen  was  then  unloaded, 
permitted  another  100  day  rest  period  and  finally  loaded  to 
a  pressure  well  in  excess  of  the  maximum  past  pressure 
experienced.   A  distinct  quasi-preconsolidation  pressure  was 
obtained.   Besides  insitu  tests  as  described  above,  consoli- 
dometer tests  were  run  on  specimens  sampled  from  the  lucite 
tube.   Good  agreement  was  obtained  provided  specimens  were 
obtained  using  a  good  piston-type  sampler  and  the  elapsed 
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time  between  unloading  and  reloading  in  the  oedometer  was 
short. 

In  the  seventh  Rankine  Lecture  of  the  British  Geotech- 
nical  Society,  Bjerrum  (1967)  introduced  his  concept  of 
delayed  consolidation  and  his  unique  relationship  between 
void  ratio,  overburden  pressure  and  time.   Figure  2.6,  taken 
from  that  paper,  illustrates  his  theory.   Line  AB  represents 
the  field  sedimentation  curve  while  BC  shows  the  change  in 
void  ratio  under  a  constant  effective  stress  due  to  delayed 
consolidation  (secondary  compression) .   If  the  soil  is  now 
loaded,  curve  CD  with  a  distinct  p   -effect  break  will  be 
obtained  and  represents  the  "instant"  (primary  consolidation 
only)  compression.   The  dotted  line  depicts  the  magnitude  of 
the  "total"  (instant  +  delayed)  compression  after  a  period 
of  time,  in  this  case  thirty  years.   This  represents  a  some- 
what different  behavior  from  that  reported  by  Leonards  and 
his  co-workers.   In  the  Bjerrum  theory  the  magnitude  of  the 
quasi-preconsolidation  pressure  is  dependent  on  the  amount 
of  secondary  compression,  while  Leonards  and  Ramiah  (1960) 
observed  a  quasi-preconsolidation  effect  even  when  secondary 
compression  was  prevented  (Figure  2.4).   Also,  in  the  void 
ratio-logarithm  pressure  plot  shown  in  Figure  2.5  (from 
Altschaeffl,  1960)  it  is  seen  that  a  much  larger  quasi- 
preconsolidation  pressure  was  observed  than  could  be  pre- 
dicted using  Bjerrum  type  curves. 


20 


a    0UD*j    piOA 


OS 

O 
tx. 


21 


2. 5  Field  Observations  of  the  Preconsolidation  Effect 
In  recent  years  several  test  embankments  have  been 
instrumented  and  the  effects  of  the  preconsolidation  pressure 
in  soft  clays  observed.   In  some  cases  it  has  been  possible 
to  take  advantage  of  this  effect  within  the  design  procedure. 
The  clay  particle  bonding  results  in  small  deflections  and 
small  excess  pore  pressures  which  dissipate  rapidly  when  the 
applied  pressure  is  less  than  some  critical  value.   When 
this  pressure  is  exceeded  deformations  are  large,  pore  pres- 
sures are  large  and  dissipation  slow. 

Such  effects  were  reported  by  H5eg  e_t  al.  (1969)  when 
two  field  load  tests  were  performed  on  a  soft  quick  Norwegian 
clay.   Piezometers  in  the  foundation  soil  showed  a  gradual 
build-up  for  low  embankment  loads  until  a  critical  stress 
level  was  attained,  after  which  a  much  higher  rate  of  in- 
crease was  observed.  Vertical  settlements  for  low  loads  in- 
creased linearly  and  were  largely  recoverable.   Beyond  the 
critical  stress  the  load  displacement  curves  started  to  bend 
over  and  the  time  dependent  settlements  accumulated  at  a 
much  higher  rate. 

Bjerrum  (1972)  as  leader  for  Session  I  -  Embankments  on 
soft  clays  at  the  Purdue  Conference  on  the  Performance  of 
Earth-Supported  Structures,  summarized  the  conference  papers 
on  preconsolidation  pressure  effects.   Papers  by  Moh  e_t  al.  , 
(1972)  and  Holtz  and  Lindshog  (1972)  demonstrate  respectively, 

behavior  when  pressures  are  less  than  p   with  pore  pressures 

cq 
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small  and  dissipation  rapid,  and  greater  than  p   when  pore 
pressures  are  large  and  dissipate  slowly.   These  effects  are 
also  evident  in  plots  of  pore  pressure  measurements  under 
low  and  high  embankments  reported  by  Elias  and  Storch  (1970)  . 
The  quasi-preconsolidation  effect  may  also  explain  results 
such  as  those  obtained  by  Rico  ejt  al.  (1969)  ,  in  which 
settlements  observed  under  three  meter  high  embankments 
were  much  smaller  than  predicted.   In  such  a  case  the  applied 
loading  probably  did  not  exceed  the  critical  pressure. 

2.6  Characteristics  of  the 
Quasi-preconsolidation  Pressure 

Compressibility :   As  has  already  been  noted  and  illustrated 
in  Figures  2.4  and  2.5,  the  existence  of  the  quasi-precon- 
solidation effect  greatly  affects  the  compressibility  of  a 
clay  after  it  has  experienced  a  period  of  time  under  con- 
stant effective  stress.  The  soil  is  capable  of  carrying 
some  additional  load  with  minimal  deformation.   Further 
increase  in  loading  will  exceed  the  critical  pressure 
rupturing  the  interparticle  bonding  and  deformations  will 
revert  to  the  virgin  compression  line. 

Coefficient  of  Consolidation:   When  a  pressure  increment  is 
applied  to  a  soil  exhibiting  a  quasi-preconsolidation  effect 
such  that  p   is  just  exceeded,  an  extremely  large  reduction 
in  the  coefficient  of  consolidation  c  is  observed.   Figure 
2.7  from  Altschaeffl  (1960)  demonstrates  this  phenomenon  for 
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two  tests  on  a  sedimented  soil.   Similar  effects  were  ob- 
served in  the  laboratory  by  Leonards  and  Ramiah  (19  60)  and 
by  HSeg,  e_t  al.  (1969)  in  the  field,  where  rates  of  pore 
pressure  dissipation  were  rapid  when  pressures  were  less 

than  p„„ .  and  slow  when  p   was  exceeded. 
rcq '  *cq 

Duration  of  Load:   Narain  et  al.  (1969)  performed  tests  on 
a  highly  plastic  organic  clay.   They  consolidated  several 
samples  to  a  pressure  of  4  tons/sq.ft.  and  allowed  the  load 
"to  remain  constant  for  different  periods  of  time  up  to  80 
days.   Continuing  the  loading  using  a  load  increment  ratio 
of  1/16  they  found  the  quasi-preconsolidation  pressure  to 
be  directly  proportional  to  the  duration  of  the  loading. 

Consolidation  Pressure :   Leonards  and  Altschaeffl  (1964) 
related  the  magnitude  of  the  quasi-preconsolidation  bonding 
to  the  initial  interparticle  spacing  and  hence  to  the  effec- 
tive overburden  pressure.   A  compilation  of  available  data 
on  the  magnitude  of  the  quasi-preconsolidation  pressure 
effect  is  shown  in  Figure  2.  8  and  shows  a  remarkable  agree- 
ment between  many  very  different  soil  deposits.   The  conclu- 
sion which  can  be  drawn  from  this  essentially  linear  plot  is 
that  for  many  soil  deposits  a  pressure  increment  of  approxi- 
mately 40  percent  of  the  effective  overburden  pressure  can 
be  applied  without  expecting  compression  to  occur  along  the 
virgin  curve.   Moh  et  al.  (1972)  report  a  higher  preconsoli- 
dation-overburden  ratio  of  1.8  for  the  soft  Bangkok  clay, 
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while  Bjerrum  (1972)  cites  a  lower  value  of  1.2-1.3  for  the 
very  soft  clays  at  Vasby,  Sweden.   Narain  e_t  al_.  (1969)  also 
found  a  linear  relationship  between  the  quasi-preconsolida- 
tion  pressure  and  the  effective  overburden  pressure  in  which 
the  slope  depended  on  the  load  increment  ratio. 

2.7  Origin  of  the  Quasi -pre consolidation  Effect 
A  variety  of  factors  has  been  proposed  to  explain  the 
physio  -  chemical  bases  for  the  quasi-preconsolidation  pres- 
sure . 

Leonards  and  Altschaeffl  (1964)  formulated  a  mechanism 
for  the  compression  of  a  clay  soil  and  for  the  formation  of 
the  quasi-preconsolidation  pressure  effect.   During  a  period 
of  time  when  a  clay  is  subjected  to  constant  applied  stresses, 
water  molecules  become  orientated  in  the  vicinity  of  the 
edge-to-face  contact  points.   Particles  slowly  displace  or 
creep  into  "the  most  efficient  arrangement  possible  from  the 
standpoint  of  bond  strength."  The  magnitude  of  this  strength 
depends  on  the  initial  interparticle  spacing  which  is  related 
to  the  effective  overburden  pressure.   The  mineral  skeleton 
can  now  sustain  pressure  increments  with  very  little  deforma- 
tion until  sliding  of  the  particles  is  again  initiated  at  the 
quasi-preconsolidation  pressure. 

Terzaghi  in  1941  postulated  the  quasi-preconsolidation 
effect  to  be  a  consequence  of  loading  rates  and  variation  of 
viscosity  of  pore  water  with  distance  from  the  particle. 


27 

Lambe  (I960)  indicates  that  particle  cementation  can 
result  from  chemical  weathering  when  ferric  oxides  which  are 
active  cementing  agents  are  present. 

Mourn  and  Rosenqvist  (1957)  suggest  another  type  of  chem- 
ical weathering  which  causes  a  change  in  the  minerology  of 
the  upper  clay  layers,  usually  the  formation  of  increased 
amounts  of  montmorillonites .   Root  holes  and  fissures  allow 
oxygen  rich  pore  water  into  the  clay  where  oxidation  of  mica 
to  montmorillonite  occurs.   Increased  plasticity  leads  to 
increased  shear  strengths. 

Bjerrum  and  Wu  (1960),  during  investigation  of  funda- 
mental shear-strength  properties  on  a  Swedish  marine  plastic 
clay,  observed  a  preconsolidation  effect  and  concluded  that 
rigid  cementation  bonds  existed  between  the  particles  in  the 
undisturbed  clay.   They  believed  the  bonds  to  be  the  result 
of  chemical  changes  over  a  geological  time,  and  possibly 
associated  with  a  period  when  the  clay  was  submerged  at 
great  depth  below  sea  level. 

Results  from  tests  performed  by  Kenny  on  a  quick  clay 
from  Labrador  indicated  strong  cementation  bonding  by  iron 
compounds,  (Bjerrum,  1967).   Samples  which  were  permeated 
with  a  chemical  which  dissolved  cementing  agents  such  as 
carbonates,  gypsum  and  iron  compounds,  showed  preconsolida- 
tion pressures  only  40  percent  of  the  magnitude  of  those 
from  untreated  specimens. 
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In  marine  clays,  the  exchangeable  ions  attached  to  the 
surface  of  the  clay  mineral  particle  are  predominantly  sodium 
cations  Na  .   Substitution  of  Calcium  Ca++,  Magnesium  Mg++, 
Aluminum  Al+++,  Ferrous  Fe++,  Ferric  Fe+++,  or  Potassium  K+ 
ions  for  the  sodium  ion  leads  to  an  increased  plasticity 
index  and  shear  strength.   The  exchange  cations  result  from 
chemical  weathering  of  felspar,  mica  or  chlorite  by  an 
acidic  pore  water. 

Thus,  it  seems  apparent  that  several  different  factors 
may  be  capable  of  producing  the  observed  quasi-preconsolida- 
tion  pressure  and  it  would  be  speculative  to  assign  a  single 
cause  for  the  formation  of  the  effect. 

2. 8   Objectives  of  Study 
From  the  review  presented,  two  important  conclusions 
may  be  drawn.   First,  the  quasi-preconsolidation  effect  is  a 
real  phenomenon  which  has  been  observed  both  in  the  labora- 
tory and  in  the  field.   It  has  a  significant  effect  on  the 
time-rate  and  magnitude  of  primary  consolidation,  causing 
normally  consolidated  deposits  to  behave  as  though  overcon- 
solidated.   Second,  it  is  noted  that  the  analyses  and  con- 
trolled experiments  conducted  in  the  past  pertain  to  one- 
dimensional  conditions  exclusively.   In  many  instances  in 
the  field  conditions  are  clearly  not  one-dimensional.   This 
suggests  a  need  to  investigate  the  quasi-preconsolidation 
phenomenon  and  its  effects  in  a  more  generalized  frame  of 
reference . 
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The  objectives  of  this  study  are,  therefore: 

1.  To  formulate  an  experimentally  derived  model 
of  the  quasi-preconsolidation  effect  under 
generalized  two-dimensional  conditions. 

2.  To  incorporate  the  quasi-preconsolidation 
effect  into  the  prediction  of  deformations 
and  stresses  during  primary  consolidation, 
for  complex  two-dimensional  boundary  and 
loading  conditions.   The  finite  element  method 
is  employed  as  the  analytical  solution  procedure 
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3.   EXPERIMENTAL  PROCEDURE  FOR  LABORATORY  TESTS 

3.1  Introduction 

In  order  to  formulate  a  two-dimensional  description  of 
the  quasi-preconsolidation  effect  a  laboratory  study  was 
performed  consisting  of  a  series  of  triaxial  compression 
•tests  on  clay  soils  in  which  a  quasi-preconsolidation  pres- 
sure (prn)  had  already  been  established.   Triaxial  test 
cq 

samples  were  cut  from  a  block  of  artificially  sedimented 
soil  and  a  quasi-preconsolidation  pressure  allowed  to  develop 
by  confining  the  samples  under  a  constant  effective  stress 
for  approximately  twenty- three  days.   This  pressure,  under 
which  the  quasi-preconsolidation  pressure  is  allowed  to 
develop,  is  denoted  herein  the  q  -  consolidation  pressure . 
The  time  during  which  it  acts  is  denoted  the  q- consolidation 
time.   Drained  compression  tests  along  a  variety  of  differ- 
ent stress  paths  were  then  performed.   Descriptions  of  the 
apparatus,  preparation  and  testing  procedures  are  given  in 
the  following  sections  and  in  Appendix  A. 

3.2  Soil  Studied 

A  commercially  available  soil,  Grundite,  processed  by 
the  Illinois  Clay  Products  Company  was  used  for  the  labora- 
tory study.   Its  clay  fraction  is  primarily  illite  (Olsen 
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and  Langfelder  1965) .   Classification  properties  are  given 
in  Figure  3.1  and  the  results  of  hydrometer  tests  in  Figure 
3.2.   The  advantage  of  using  such  a  commercial  clay  is  that 
large  quantities  with  minimal  variation  of  the  material 
properties  can  be  obtained. 

3 .3  Sedimentation  Samples 
The  clay  samples  to  be  ultimately  tested  under  triaxial 
conditions  were  initially  prepared  by  a  one-dimensional 
sedimentation  process  in  3  1/2  inch  diameter  lucite  tubes. 
The  loading  frame,  in  which  very  slow  rates  of  sedimentation 
can  be  achieved,  and  the  sedimentation  equipment  are 
described  in  Appendix  A.l.   All  samples  were  prepared  in  an 
identical  manner,  using  the  same  proportions  of  grundite, 
water,  dispersant  and  flocculant.   Mixing  and  deairing 
periods  were  also  constant.   Sedimentation  samples  were 
initially  loaded  in  the  frame  at  a  very  slow  rate  which  was 
then  increased  such  that  a  final  pressure  of  11.45  psi  was 
reached  in  approximately  three  weeks.  Tests  on  the  homo- 
geneity of  the  samples  were  made  and  are  described  along 
with  a  detailed  description  of  the  preparation  procedure  in 
Appendix  A. 2. 

3.4  Triaxial  Samples 
Triaxial  test  specimens  of  standard  dimensions  were  cut 
from  the  sedimented  soil  described  in  the  previous  section. 
Details  of  the  procedure  are  given  below.   It  was  then 
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GRUNDITE     SOIL 
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desired  to  isotropically  load  the  sample  and  maintain  a 
constant  cell  confining  pressure  for  a  period  of  approximately 
three  weeks.   A  set-up  whereby  eight  triaxial  cells  could  be 
simultaneously  loaded  was  designed.   The  confining  pressure 
was  supplied  via  a  laboratory  air  source  with  a  back  up 
nitrogen  tank  for  use  should  the  primary  source  fail  to 
maintain  the  desired  pressure.   This  triaxial  cell  set-up 
is  described  in  detail  in  Appendix  A. 3.1.   Details  of  the 
triaxial  cells,  the  constant  strain  rate  testing  machine 
and  the  recorders  used  are  given  in  Appendices  A. 3. 2  and 
A. 3. 3. 

The  laboratory  sedimented  sample  was  extruded  from  the 
lucite  tube  by  removing  the  base  and  porous  stone  and  push- 
ing the  sample  out  with  the  piston.   The  sample  was  approxi- 
mately 4  inches  high,  with  a  unit  weight  of  112  pcf  and  a 
void  ratio  of  1.57.   The  sample  was  trimmed  using  the  soil 
lathe  shown  in  Figure  3-3  to  a  right  circular  cylinder,  1.4 
inches  in  diameter,  and  2.9  inches  high.   Slotted  filter 
paper  was  placed  around  the  sample  to  facilitate  drainage, 
and  deaired  saturated  porous  stones  were  placed  on  each  end. 
The  sample  with  top  platen  was  mounted  on  the  cell  pedestal 
and  enclosed  with  two  rubber  membranes,  held  in  place  by 
O-Rings  or  strips  of  Dental  Dam.   Drainage  lines  were  all 
carefully  filled  with  water  taken  from  the  sedimentation 
tube  with  the  same  chemical  composition  as  the  sample  pore 
fluid.   The  triaxial  cell  was  screwed  into  place  and  filled 
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FIGURE     3-3    TRIAXIAL    SPECIMEN     PREPARATION     EQUIPMENT. 
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with  distilled  water.   A  1/2  inch  layer  of  S.A.E.  40  oil  was 
placed  at  the  top  of  the  cell  to  insure  piston  lubrication 
during  the  long  loading  period. 

The  triaxial  cell  was  then  connected  into  the  system 
and  a  burette  attached  to  the  drainage  lines.   The  confining 
pressure  was  increased  over  a  period  of  two  days  to  the 
desired  value  and  the  sample  allowed  to  remain  under  this 
constant  pressure  for  approximately  twenty  three  days. 
Burette  levels  were  read  and  recorded  regularly. 

3.5  Definitions 

Prior  research  on  the  quasi-preconsolidation  effect  in 
clay  soils  has  been  carried  out  under  conditions  of  one 
dimensional  loading  and  drainage.   In  such  a  test  only  one 
stress  parameter,  the  vertical  pressure  p,  and  one  strain 
parameter,  the  vertical  strain  e  or  the  void  ratio  e,  are 
required  for  complete  recording  of  the  test  behavior, 
usually  in  the  form  of  an  e  -  log  p  plot  (Leonards  and 
Altschaeffl   1964)  . 

In  this  work  the  effect  has  been  studied  using  the  more 
versatile  triaxial  test  apparatus,  in  which  both  the  vertical 
and  radial  stresses  can  be  varied.   It  is  assumed  that  these 
are  respectively  equal  to  the  major  and  minor  principle 
stresses  and  that  the  tangential  stress  is  always  the  inter- 
mediate principle  stress  and  is  equal  in  magnitude  to  the 
radial  stress  (Perloff  and  Pombo   1969)  .   The  axial  stress 
is  equal  to  the  total  vertical  load  divided  by  the  sample 
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area  while  the  radial  stress  is  equal  to  the  confining  fluid 
pressure  A  diagram  illustrating  the  triaxial  test  stress 
state  is  given  in  Figure  3.4. 

Stress  parameters  used  for  the  representation  of  the 
triaxial  test  results  were  based  on  those  recommended  by 
Roscoe  and  Burland  (1968)  and  are  derived  below. 

A  stress  increment  is  defined  as  the  difference  between 
the  final  and  the  initial  stress  levels 

i.e.  La.    =   a.  -  o.  j  =  1,  2,  3       (3.1) 

■*    J  final   ■'initial 


or   Aa.  =  a.      -  o. 

Octahedral  Normal  Effective  Stress  o\  «. 

oct 

aoct  =  (0J  +  o 2  +  o,)/3  (3.2) 

°OCt  ~  aOCtj:.    ,   "  aOCt.„.  +  .  , 

final       initial 


■  Oi£  +  a2£  +  o3f)/3  -  (a^  +  a2i  +  o3;L)/3 

=  (AcTj  +  Aa2  +  Aa3)/3  (3.3) 

For  triaxial  test  results  this  reduces  to 

Aaoct  =  (Aa,  +  2  Ao"3)/3  (3.4) 
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Octahedral  Shear  Stress 

Toct  =  '(oTj  -  02)2  +  (c2  -  03)2  +  (o,  -  a,)2  (3.5) 

At     =  t  -  T 

OCt     OCtj--   n     oct.  ...  i 
final       initial 

■  /(alf  -  o*2£)2  +  (a2£  -  o3f)2  +  (o"3f  -  olf)2 

-  /(ali  -  a2i)2  +  (a2i  -  03j_)2  +  (a3±   -   o1±)2    (3.6) 
For  triaxial  test  results  this  reduces  to 


A 


Toct  =  3  (A°i  "  Aa3^  C3-7) 


A  useful  parameter  for  designating  certain  stress  paths  is 
the  octahedral  effective  stress  ratio. 

oct  /.-  „•> 

n  »  — (3.8) 

oct 
As  with  the  stress  state,  three  principle  natural 
strains  can  be  defined.   In  compression  tests  the  axial  and 
radial  strains  are  respectively  the  major  and  the  minor 
principle  strains,  with  the  assumption  that  ex  >  e2  =  e3. 
The  two  strain  increment  parameters  chosen  were  the  com- 
pressive volumetric  strain  increment  5v  and  the  deviatoric 
(shear  distortional)  strain  increment  fie,  where 

6v  =  -  ^  =  -  jA^-  =  5£l  +  fie2  ♦  6e3         (3.9) 
fie  =  6ex  -  6v/3  (3.10) 
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where  AV  and  Ae  are  the  changes  in  volume  and  void  ratio. 
For  the  triaxial  test  these  reduce  to 

<5v  =  6e2  +  2  5e3  (3.11) 

and  6e  =  I  (6e1    -   6e3)  (3.12) 

3.6  Triaxial  Testing 


Triaxial  tests  along  four  different  types  of  stress 
path  were  carried  out.   These  paths  are  illustrated  in 
Figure  3.5  and  are  designated  as  follows 

Type  1   Constant  chamber  pressure  tests 

Type  2   Isotropic  loading  tests 

Type  3  Constant  o      .    tests 

Type  4   Constant  n  tests 
The  constant  chamber  pressure  tests,  Type  1,  were  run 
in  the  Engineering  Laboratory  Equipment  (E.L.E.)  frame 
described  in  Appendix  A. 3. 3.   In  the  special  case  where  the 
test  chamber  pressure  was  equal  to  q  ,  the  confining  pressure 
under  which  the  quasi-preconsolidation  effect  was  allowed  to 
develop,  an  extension  piece  of  tubing  was  inserted  between 
the  cell  and  the  triaxial  set-up  and  the  cell  placed  in  the 
E.L.E.  frame.   When  the  test  confining  pressure  was  differ- 
ent from  q  ,  the  cell  was  connected  independently  to  one  of 
the  air  sources  and  the  chamber  pressure  changed  slowly  in 
increments  over  a  period  of  several  hours  or  days,  to  the 
desired  test  pressure.   The  load-cell  was  carefully  calibrated 
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FIGURE     3-5      TRIAXIAL   TEST  STRESS    PATHS 


42 

using  dead  weight  and  the  sample  volume  change  burette  also 
calibrated.   This  simple  apparatus  shown  in  Figure  3.6  con- 
sisted of  a  2  cc  burette  connected  to  the  sample  drainage 
lines.   A  pressure  transducer  measured  the  weight  of  the 
column  of  water  and  the  output  was  recorded  on  the  same  chart 
as  the  load  cell  readings.   Most  of  the  tests  were  run  at  a 
constant  rate  of  4.0  x  10   inches  per  minute,  which  is 
approximately  one  percent  strain  in  twelve  hours.   At  the 
end  of  a  test  water  contents  at  six  different  levels  within 
the  sample  were  obtained. 

A  computer  program  was  written  to  reduce  the  triaxial 
test  data,  and  is  described  in  detail  in  Appendix  C  where 
input  format,  flow  chart,  program  listing  and  an  example 
problem  are  given. 

The  Isotropic  loading  tests  (Type  2)  were  carried  out 
in  the  adapted  Geonor  cell  described  in  Appendix  A. 3. 2  and 
illustrated  in  Figure  A. 15.   After  the  q  confinement  period 
the  Geonor  cell  was  independently  connected  to  one  of  the 
sources  and  the  confining  pressure  increased  incrementally 
at  regular  intervals.   Burette  levels  were  noted  at  the  end 
of  each  loading  period  and  output  from  the  LVDT  was  continu- 
ously recorded  on  the  Varian  machine. 

The  constant  a   t  and  constant  n  tests  (Types  3  and  4) 
were  run  using  standard  Geonor  cells  in  a  special  incre- 
mental loading  frame.   The  set-up  is  shown  in  Figure  3.7. 
The  cell  confining  pressure  was  supplied,  as  in  all  the 
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FIGURE    3-6     BURETTE     WITH     PRESSURE     TRANSDUCER 
ATTACHMENT . 
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FIGURE     3-7    INCREMENTAL    LOADING    FRAME. 
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tests,  by  compressed  air,  and  controlled  by  means  of  a 
regulator  with  attached  dial  gauge.   The  axial  stress  was 
applied  incrementally  by  adding  weights  to  an  initially 
balanced  hanger  bearing  on  the  cell  piston.   For  the  constant 

5L,.*  stress  path  it  is  necessary  that 

oct       r  * 

-    .  Ao,  +  2   Aa3  =  0 
oct       3 

or 

Aa3  =  -  Aa1/2  (3.13) 

This  test  was  conducted  by  increasing  the  axial  stress  by 
some  increment  and  decreasing  the  confining  pressure  by  half 
as  much.   In  Type  4  tests  similar  computations  were  carried 
out  for  each  n  value.   For  example 


oct 


.  .       J-j  (Aox  -  Aa3)  =  j   (Aax  +  2  Aa3) 

or  Aa3  =  0.121  Aal 

In  this  case  the  confining  pressure  was  increased  by  an 
amount  equal  to  0.121  times  the  axial  stress  increment. 
A  Type  3  test  is,  in  fact,  a  special  case  of  Type  4  with 
ti  =  °°.   Since  pistons  of  smaller  cross  sectional  area  than 
the  sample  were  used  the  load  L  required  to  give  a  desired 
vertical  stress  a     is  given  by 
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L  =  a  A^  -  a   (Ao  -  A  ) 
V  s    c  v  s    py 


where   a   is  the  cell  confining  pressure 
A  is  the  sample  area 
A  is  the  piston  area. 
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4.   RESULTS  AND  DISCUSSION  OF  QUASI  TESTS 

4. 1   Introduction 

Twenty-seven  triaxial  test  specimens  were  prepared  and 
allowed  to  consolidate  under  constant  confining  pressures  q_ 
as  described  in  Section  3.4.   Ten  samples  were  confined  under 
a  constant  pressure  of  16.7  psi,  eight  under  24.0  psi  and  the 
remaining  nine  under  30.0  psi.   This  period  of  isotropic  con- 
solidation, called  the  q-consolidation  time,  ranged  from  20 
to  27  days  in  all  cases  except  test  Number  21,  which  was 
allowed  to  remain  44  days  prior  to  testing. 

All  the  grundite  sedimentation  samples  were  prepared  in 
the  same  manner  and  all  but  Test  Number  1  loaded  to  50  kg.  in 
the  sedimentation  frame.   The  resulting  specimens  had  a  void 
ratio  range  of  1.52  to  1.62.   In  the  case  of  Test  1  the  sedi- 
mentation sample  was  loaded  to  30  kg.,  with  a  resulting 
higher  void  ratio  of  1.88. 

Three  types  of  triaxial  cells  were  used  and  are  described 
in  Section  A. 3. 2.   Twenty-two  tests  were  performed  in  Geonor 
cells,  four  in  the  Wykeham-Farrance  and  one  in  the  Testlab 
cell.   Five  tests  were  carried  out  with  single  drainage  from 
the  base  of  the  sample,  while  the  remaining  twenty- two  had 
both  top  and  bottom  drainage  facilities.   Test  details  are 
given  in  Table  4.1. 


TABLE  4.1 
q-Consolidation  Pressures  and  Times 
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Test 

1 

2 

3 

4 

5 

6 

7 

8 

9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 


Cell 

Double 

or 

Single 

Drainage 

Void 
Ratio 

qc 

qc  Time 
(days) 

Testlab 

Single 

1.88 

16.7 

21 

Wykeham 

Double 

1.53 

16.7 

21 

Geonor 

Double 

1.55 

16.7 

24 

Geonor 

Double 

1.54 

16.7 

27 

Geonor 

Double 

1.61 

16.7 

25 

Geonor 

Double 

1.60 

16.7 

26 

Geonor 

Double 

1.62 

16.7 

20 

Geonor 

Single 

1.52 

16.7 

22 

Geonor 

Double 

1.61 

16.7 

22 

Geonor 

Double 

1.57 

16.7 

20 

Geonor 

Double 

1.53 

24.0 

25 

Geonor 

Double 

1.57 

24.0 

25 

Wykeham 

Double 

1.53 

24.0 

24 

Geonor 

Double 

1.55 

24.0 

21 

Geonor 

Single 

1.56 

24.0 

22 

Geonor 

Double 

1.55 

24.0 

26 

Geonor 

Double 

1.56 

24.0 

22 

Geonor 

Double 

1.61 

24.0 

25 

Geonor 

Double 

1.54 

30.0 

21 

Geonor 

Double 

1.54 

30.0 

27 

Geonor 

Single 

1.61 

30.0 

44 

Geonor 

Double 

1.60 

30.0 

21 

Wykeham 

Double 

1.57 

30.0 

22 

Geonor 

Double 

1.55 

30.0 

25 

Wykeham 

Single 

1.60 

30.0 

24 

Geonor 

Double 

1.62 

30.0 

22 

Geonor 

Double 

1.56 

30.0 

21 
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The  test  results  are  summarized  and  interpreted  in  the 
following  sections. 

4 . 2   The  Triaxial  Test  Types 
Four  different  types  of  triaxial  tests  described  in 
Section  3.6  were  performed. 

In  summary,  Type  1  were  constant  strain  rate  tests  with 
constant  confining  pressures.   Type  2  were  isotropic  tests 
using  the  adapted  geonor  cell  with  LVDT  attachment  in  which 
the  confining  pressure  was  increased  incrementally.   Type  3 
were  incremental  loading  tests  in  which  the  confining  pres- 
sure was  decreased  and  the  axial  pressure  increased  so  that 
a  constant  octahedral  normal  stress  was  maintained.   Type  4 
were  also  incremental  loading  tests  in  which  the  axial  and 
confining  stresses  were  altered  so  as  to  maintain  a  constant 
octahedral  shear  stress,  normal  stress  ratio  n- 

Table  4-2  lists  the  test  numbers,  types  of  tests  per- 
formed, rates  of  loading,  and  qc.,  the  octahedral  normal 
stress  prior  to  increasing  the  shear  stress. 

4.3   The_  q-Yield  Surface 
The  one-dimensional  consolidation  testing  in  which 
quasi-preconsolidation  effects  were  observed  has  been  dis- 
cussed in  Chapter  2  and  typical  e-log  p  plots  of  the  behavior 
are  shown  in  Figures  2.4  and  2.5.   The  salient  feature  of  the 
testing  was  the  apparent  formation  of  interparticle  bonding 
during  a  rest  or  confinement  period  at  constant  effective 
stress,  which  enabled  the  soil  to  carry  an  additional  loading 
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TABLE  4.2 
Test  Types  and  Loading  Rates 


Test 

Type 

<*C 

q'i 

1 

1 

16.7 

16.7 

2 

1 

16.7 

16.7 

3 

1 

16.7 

16.7 

4 

2 

16.7 

- 

5 

1 

16.7 

16.7 

6 

2 

16.7 

- 

7 

2 

16.7 

- 

8 

1 

16.7 

16.7 

9 

1 

16.7 

18.2 

10 

1 

16.7 

19.7 

11 

1 

24.0 

30.0 

12 

2 

24.0 

- 

13 

1 

24.0 

24.0 

14 

3 

24.0 

24.0 

15 

3 

24.0 

25.5 

16 

4 

24.0 

24.0 

17 

1 

24.0 

25.0 

18 

3 

24.0 

30.0 

19 

1 

30.0 

30.0 

20 

1 

30.0 

35.0 

21 

1 

30.0 

30.0 

22 

2 

30.0 

- 

23 

1 

30.0 

30.0 

24 

1 

30.0 

40.0 

25 

1 

30.0 

26.0 

26 

4,3 

30.0 

30.0 

27 

4 

30.0 

30.0 

Loading  Rate 

5.0  x  10_lf  inches/minute 
1.5  x  10"1*  inches/minute 
4.0  x  10" 5  inches/minute 
0.5  psi/2  hours 
1.5  x  10"1*  inches/minute 
0.5  psi/3  hours 
1.0  psi/hour 

1.5  x  io "^  inches/minute 
4.0  x  10"1*  inches/minute 
1.5  x  10"5  inches/minute 
4.0  x  10"5  inches/minute 
1.0  psi/12  hours 
4.0  x  10" 5  inches/minute 
increase  Ao,  1  psi/24  hours 
increase  Aoj  1  psi/24  hours 
increase  Ao^  1  psi/24  hours 
4.0  x  10" 5  inches/minute 
increase  AOj  1  psi/24  hours 
4.0  x  io" 5  inches/minute 
4.0  x  io-5  inches/minute 
4.0  x  10"5  inches/minute 
1.0  psi/12  hours 
4.0  x  io-5  inches/minute 
4.0  x  10" 5  inches/minute 
4.0  x  10" 5  inches/minute 
decrease  Aa3  5 psi/24  hours 
increase  AOj  4  psi/24  hours 
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with  only  small  deformation.   If  this  loading  were  increased 
beyond  some  critical  pressure,  the  quasi-preconsolidation 
pressure,  the  bonds  were  ruptured  and  the  e-log  p  plot  reverted 
to  the  virgin  curve.   This  was  shown  schematically  in  Figure 
2.2.   On  Figure  2.8  are  plotted  the  results  of  the  three  iso- 
tropic loading  tests  for  different  q-confining  pressures. 
Although  not  directly  comparable  because  of  the  isotropic 
q-conf inement  pressure  the  results  do  compare  well  with  this 
field  data. 

Figures  B.4  to  B.41  for  the  triaxial  tests  performed  in 
this  study  indicate  a  similar  type  of  behavior.   The  origin 
in  each  plot  represents  the  octahedral  normal  stress  qc .  in 
the  triaxial  cell  prior  to  increasing  the  shear  stress.   The 
axes  denote  stress  increments  from  this  "initial"  state.   For 
all  the  tests  performed  stress-strain  curves  are  initially 
steep,  i.e.,  the  soil  exhibits  low  compressibility,  then  change 
fairly  abruptly  to  a  more  shallow  (more  compressible)  slope. 
All  plots  have  the  appearance  of  two  straight  lines  with  a 
transition  curve,  until  large  strains  when  failure  ensues. 

For  each  test  the  intersection  of  the  two  straight  line 
portions  was  defined  as  the  q-yield  stress  point.   The  vertical 
and  lateral  stresses  at  this  point  are  known  and  from  these  the 
octahedral  normal  (Ao7-,,*.  )  and  shear  (At__+  )  stresses  can  be 
calculated.   These  are  listed  in  Table  4.3,  and  are  plotted 
along  with  Test  Numbers  and  Stress-paths  in  Figures  4.1,  4.2, 
or  4.3,  each  figure  corresponding  to  a  specific  q-consolidation 
pressure.   Each  plot  depicts  a  curve  which  will  be  called  the 
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TABLE  4.3 

q-Yield  Stress  Results 

Test 
Number 

i0°ctq 

AToct 

q 

1 

2.75 

3.89 

2 

2.80 

3.96 

3 

2.50 

3.54 

4 

6.10 

0.00 

5 

2.90 

4.10 

6 

6.50 

0.00 

7 

6.30 

0.00 

8 

2.73 

3.86 

9 

4.03 

3.58 

10 

5.20 

3.11 

11 

8.20 

3.11 

12 

9.00 

0.00 

13 

3.75 

5.30 

14 

0.00 

5.15 

15 

1.50 

5.40 

16 

-4.80 

1.20 

17 

4.47 

4.90 

18 

6.00 

4.60 

19 

5.23 

7.40 

20 

8.89 

5.50 

21 

4.66 

6.59 

22 

11.50 

0.00 

23 

5.10 

2,20 

24 

11.50 

2.12 

25 

0.81 

6.80 

26 

2.88 

7,20 

'27 

11.30 

2.83  - 
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q-yield  curve  for  that  particular  q-consolidation  pressure. 
A  Stepwise  Regression  Program  BMD2R,  Purdue  University 
Computer  Science  Center,  was  used  to  calculate  the  best-fit 
polynomials  for  the  three  sets  of  quasi  data  points.  The 
resulting  curves  were  ellipses  displaced  along  the  abscissa 
and  having  equations  of  the  form 

y2  =  Cl    +  C2  x  +  C3  x2  (4-1) 

The  values  of  the  constants  are  given  in  Table  4.4.   The 
"Best-fit"  curves  have  been  plotted  with  the  data  points  in 
Figures  4.1,  4.2  and  4.3,  and  are  drawn  together  in  Figure 
4.4.   It  was  then  desired  to  formulate  a  single  equation 
which  would  represent  the  q-yield  curve  for  any  q-consolida- 
tion pressure.   The  method  is  illustrated  schematically  in 
Figure  4.5.   In  Figure  4.4  straight  lines  radiating  from  the 
origin  at  angles  a  (0,  5,  10,  20,  30,  40,  54.7  and  90°)  were 
constructed.   The  ordinates  (y  )  of  the  points  of  intersec- 
tion of  these  a-lines  and  the  three  q-yield  curves  were 
noted.   For  the  a  =  0  line  where  the  intercept  ordinates  are 
zero  the  abscissae  values  (x  )  were  recorded  and  used  instead, 
In  Figure  4.6  these  intercept  values  were  plotted  on  their 
respective  q-consolidation  pressure  lines,  16.7,  24.0  and 
30.0  psi.   Straight  lines  known  as  g-lines  were  then  drawn 
through  the  points  representing  a  particular  a.   A  least 
squares  analysis  was  used  to  determine  the  best  straight 
lines,  and  their  slopes  and  intercept  values  with  a  20  psi 
consolidation  pressure  line  are  given  in  Table  4.5.   Since 
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TABLE  4.4 
Ellipse  Constants  for  the  q-Yield  Curves 


q- Consolidation 

Pressure 

psi 

c, 

c2 

c3 

16.7 

12.424 

3.310 

-0.823 

24.0 

26.477 

2.489 

-0.583 

30.0 

42.039 

4.746 

-0.703 

TABLE  4.5 
Least-squares  Parameters  for  8-lines 


Intercept  with 

a-line  angle 

Slope  of  0-line 

20 

psi  line 

0° 

0.3885 

7.769 

5° 

0.0338 

0.676 

10° 

0.0661 

1.322 

20° 

0.1214 

2.428 

30° 

0.1620 

3.240 

40° 

0.1892 

3.783 

55° 

0.2280 

4.560 

90° 

0.2149 

4.297 
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these  B-lines  all  intersect  at  a  zero  consolidation  pressure 
point  the  q-yield  curves  represent  sections  of  an  elliptic 
cone,  shown  sketched  in  Figure  4.7,  and  called  the  q-yield 
surface.   The  q-yield  curve  for  a  grundite  specimen  prepared 
under  any  q-consolidation  pressure  can  now  be  found  by  inter- 
polation from  the  known  curves.   The  intercepts  of  the  3-lines 
and  a  20  psi  consolidation  line  were  observed  in  Figure  4.6 
and  noted  in  Table  4.5.   These  points  were  used  to  calculate 
the  20  psi  interpolation  q-yield  curve 

y2  =  18.903  +  1.617  x  -  0.526  x2         (4-2) 

A  study  of  the  form  of  this  equation  with  relation  to  an 
elliptic  cone  leads  to  the  general  equation 

y2  =  0.032818  q2  +  0.067384  qx  -  0.52570  x2     (4-3) 

where 

q  is  the  q-consolidation  pressure  (psi) 

x  is  the  q-octahedral  normal  stress  increment  (psi) 

y  is  the  q-octahedral  shear  stress  increment  (psi) 

Use  will  be  made  of  this  equation  later  in  the  Finite 
Element  formulation  in  Chapter  5. 
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4. 4   Discussion 

The  aim  of  the  laboratory  experimentation  was  to  formu- 
late a  simple  soil  model  which  would  take  into  account,  in 
a  more  general  manner,  the  quasi-preconsolidation  effect 
observed  by  Leonards  and  Ramiah  (1960)  for  one-dimensional 
consolidation  tests.   This  was  accomplished  by  running  tri- 
axial  tests  on  clay  samples  which  had  been  prepared  in  a 
manner  conducive  to  the  formation  of  quasi-preconsolidation 
bonding.   The  samples  were  prepared  and  the  tests  run  as 
described  in  Chapter  3  and  Appendix  A.   The  resulting  stress- 
strain  curves  showed  an  initially  steep  slope  followed  by  a 
fairly  sharp  transition  to  an  essentially  linear  curve  of 
much  lower  slope.   The  intersection  of  the  two  straight  line 
portions  of  the  curve  was  termed  the  q-stress  point,  and 
represents  the  stress  at  which  the  temporary  bonding  which 
had  developed  during  the  consolidation  period  was  assumed 
to  rupture.   Bond  breaking  probably  occurred  over  some  range 
of  stresses,  explaining  the  transition  curve.   Conlon  (1966) 
observed  similar  shaped  curves  for  cemented  Canadian  marine 
and  estuarine  clays.   He  explained  the  shape  as  the  sum  of 
two  effects,  namely,  the  normal  soil  shearing  resistance  and 
the  interparticle  bonding  strength.   This  is  shown  schematic- 
ally in  Figure  4.8. 

The  q-stress  points  for  different  tests  and  different 
q-consolidation  pressures  were  plotted  to  form  q-yield  curves 
and  the  q-yield  surface  as  described  in  Section  4.3.   Of  the 
twenty-seven  tests  there  were  19  different  stress  paths;  the 
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FIGURE     4-8 


SCHEMATIC  DEVELOPMENT  OF  THE  SHAPE 
OF  THE  QUASI  STRESS- STRAIN  CURVES 
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remaining  8  tests  being  duplicate  tests  mostly  of  Type  1  with 
confining  pressures  equal  to  their  q-consolidation  pressure. 
With  such  a  limited  number  of  points  spread  over  an  extensive 
three-dimensional  surface  it  is  possible  only  to  state  that 
there  appears  to  be  a  unique  surface,  independent  of  stress 
path,  which  divides  this  stress  space  into  two  regions,  A  and 
B.   Space  B  within  the  surface  and  A  outside  the  surface  repre- 
sent soil  behavior  Before  and  After  rupture  of  the  quasi-pre- 
consolidation  bonds.   Since  the  stress-strain  plots  showed 
■essentially  bilinear  behavior  a  simple  soil  model  may  be 
defined  giving  the  soil  a  large  Modulus  value  Eg  if  its  stress 
path  lies  within  Region  B.   When  a  stress  path  crosses  the 
surface  the  bonds  are  assumed  to  rupture  and  a  lower  Modulus 
E^  would  be  appropriate.   The  model  proposed  is  entirely  elastic 
and  no  consideration  is  made  of  ultimate  yield.   The  inherent 
assumption  is  therefore  made  that  loads  do  not  exceed  the 
elastic  limit  and  no  plastic  behavior  is  permissible.   The 
effects  of  the  q-yield  surface  and  the  bilinear  constitutive 
law  during  the  consolidation  process  will  be  demonstrated  in 
Chapter  5. 

During  the  sample  preparation  and  testing  little  attempt 
was  made  to  consider  the  effects  of  all  the  parameters  which 
might  change  the  magnitudes  of  the  q-stress  points.   In  fact 
every  effort  was  made  to  standardize  the  complete  testing 
procedure.   This  was  not  always  possible  especially  for  the 
q-consolidation  time  since  actual  testing  took  varying 
amounts  of  time  and  the  equipment  was  not  always  vacant  for 
a  particular  sample  at  a  particular  time.   Table  4.1  shows 
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that  the  consolidation  periods  varied  from  20  to  27  days. 
In  order  to  ascertain  if  this  time  might  be  a  critical  param- 
eter one  test,  Number  21,  was  permitted  to  remain  under  a 
consolidation  pressure  of  30  psi  for  44  days,  a  period  almost 
twice  as  long  as  the  average  value  of  23.9  days.   It  might 
have  been  expected  that  since  more  time  was  available  for 
inter  particle  movement  and  bonding,  a  higher  q-stress  point 
would  result,  as  was  observed  by  Narain  et  a_l.  (1969)  for 
the  one-dimensional  case.   This  was  not  so;  in  fact  test 
Number  21  had  the  lowest  q-stress  point  of  the  three  samples 
tested  along  that  particular  stress  path.   Since  this  test 
did  not  show  any  increase  in  q-stress,  let  alone  the  almost 
2  to  1  increase  obtained  by  Narain  in  one-dimensional  testing, 
it  may  be  assumed  that  the  effects  of  the  length  of  q-consoli- 
dation  with  a  maximum  of  4  days  variation  from  the  mean  may 
be  neglected. 

Triaxial  Test  Number  1  was  loaded  in  the  sedimentation 
tube  to  a  final  pressure  40  percent  less  than  all  other  tests 
and  consequently  had  a  higher  initial  void  ratio.   However 
during  the  q-consolidation  period  it  should  attain  the  same 
void  ratio  as  the  other  samples  assuming  a  unique  e-log  p 
relationship  for  the  grundite,  and  this  was  apparently  borne 
out  since  the  quasi  stress  for  Test  1  was  approximately  the 
mean  of  the  five  similar  tests. 

Table  4.1  shows  that  while  most  of  the  triaxial  speci- 
mens had  drainage  from  both  top  and  bottom,  six  tests  had 
bottom  or  pedestal  drainage  only.   Table  4.3  shows  that  most 
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of  the  tests  performed  in  the  Testing  Apparatus  were  run  at 
a  rate  of  4.0  *  10" 5  inches/minute  while  a  few  were  run  at 
1.5  x  10""  and  one  at  5.0  x  10"14  inches/minute.   These  two 
factors  can  be  lumped  together  as  a  measure  of  the  ease  with 
which  pore  water  may  flow  out  of  the  sample  so  as  to  maintain 
a  zero  or  low  excess  pore  pressure  and  a  pure  effective 
stress  system.   Triaxial  Test  Number  3  had  double  drainage 
and  a  loading  rate  of  4.0  x  10" 5  inches/minute  while  Test 
Number  8  had  single  drainage  and  the  faster  loading  rate  of 
1.5  x  10"1*  inches/minute.   Other  factors  were  essentially 
similar.   Test  8  resulted  in  a  q-octahedral  shear  stress 
about  8  percent  higher  than  in  Test  3.   This  increase  is  in 
the  direction  reported  by  Crawford  (1964,  1965)  during 
studies  on  the  effect  of  loading  rate  on  the  preconsolidation 
pressure.   A  recommended  rate  of  loading  can  be  calculated  by 
a  method  formulated  by  Gibson  and  Henkel  (1954) .   They  found 
that  the  average  degree  of  dissipation  at  failure  Ur  may  be 
expressed  in  the  form 

trf  -  i  -  ^  (4-4) 

where  h  =  half  the  height  of  the  sample 

c  =  coefficient  of  consolidation 

tr  =   time  to  failure 
and    5  =  a  factor  depending  on  drainage  conditions. 

For  a  degree  of  dissipation  of  951,  the  time  to  failure  is 
given  by 
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20  h2 
t£  "  TcT"  (4-5) 


For  a  test  with  drainage  from  both  ends  and  filter-paper 
drains  Bishop  and  Henkel  (1962)  recommend  a  value  of  £  =  40.4. 
The  value  of  c  was  obtained  from  a  consolidation  test  per- 
formed on  Grundite.   From  a  dial  reading  —  time  plot  for  a  load 
increment  beyond  the  quasi  preconsolidation  yield  point,  t 

was  found 

_   .  0.197  H2 

v  "   t77~ 

(4-6) 
■  6.94  x  10" H  in?/min. 


The  time  to  failure 

t,  =  2(?  X  -1;42 r-  =  1398  min.      (4-7) 

f   40.4(6.94  x  10"H) 

If  failure  occurs  at  15  percent  strain  then  the  strain  or 
loading  rate  is  given  by 

0.028  x 15 
erate  =  1398  (4-8) 

=  3.00  x  10""  in. /min. 

With  the  exception  of  Test  Number  1  this  recommended  rate  is 
faster  than  the  two  loading  rates  (1.5  x  10"1*  and  4.0  x  10"5 
in. /min.)  used  for  the  triaxial  testing. 
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5.   THE  FINITE  ELEMENT  SOLUTION  PROCEDURE 

5. 1   Introduction 

The  finite  element  program  into  which  the  quasi-precon- 
solidation  soil  model  has  been  introduced  is  an  adaptation 
of  a  consolidation  program  CONSOL  written  at  the  Massachusetts 
Institute  of  Technology,  (Christian  and  Boehmer,  1970; 
Christian  ejt  a_l .  ,  1972),  and  based  on  programs  developed 
originally  by  Wilson  (1963) . 

Single  phase  finite  element  programs  have  been  widely 
reported  in  the  Soil  Mechanics  literature  and  formulation  of 
the  equations  is  readily  available  (Perloff,  1969).   However, 
since  the  consolidation  process  is  a  two  phase  problem  (soil 
skeleton  and  pore  water) ,  a  short  description  of  the  theory 
and  equation  derivation  is  given  in  the  following  two  sec- 
tions.  Section  5.4  briefly  describes  the  original  CONSOL 
program  setup,  changes  incorporated,  and  the  introduction  of 
the  quasi-preconsolidation  model  formulated  from  the  labora- 
tory experiments.   In  Section  5.5  results  from  example 
programs  illustrate  the  effects  of  the  model  on  pore  pressures, 
displacements  and  stresses  during  the  primary  consolidation 
process.   In  the  final  section  qualitative  comparison  is  made 
between  analytical  results  and  recently  published  field 
studies. 
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5.2   Finite-element  Equations  for  a  Two-phase  Material 
The  principle  of  virtual  displacements  states  that  for 
a  conservative  system  under  equilibrium  the  virtual  work  dW, 
due  to  any  virtual  displacement  from  a  compatible  state  is 
equal  to  the  virtual  strain  energy  dU.. 


dW  =  dU 


(5-1) 


where   dW  is  the  first  order  work  done  by  the  external  forces 
dU^  is  the  first  order  work  done  by  the  internal 
stresses . 
The  internal  work  is  equal  to 


dU 


6{e}   {a}  d(volume) 


where  the  strain  vector  is 
T 


{£}i  =  [6ex  6ey  6^xy_ 


and  the  stress  vector 
[o]   = 


xy 


The  external  work  is  equal  to 


6W 


(5-2) 


(5-3) 


(5-4) 


=  f  6{u}   {b}  dv  +  f  6{u}   {p}  ds  +  6{u  }e'T  {F}*    (5-5) 
Jv  '  s  n        c 


Equating  the  external  work  and  the  internal  strain  energy  and 
substituting  in  the  effective  stress  principle 
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{a}   =  {a}  +  {S}  (5-6) 

gives 


6{e}T  {a}  dv  + 


6{e}T  {S}  dv  = 


6{u}T  {b}  dv 


[  6{u}T  {p}  ds  +  6{un}e'T  {T<}ec  (5-7) 


+ 
s 

where 

{a}   the  effective  stress  associated  with  the  soil 

skeleton 

{S}   the  excess  pore  pressure  associated  with  the  pore 

fluid  phase 

{b}   the  external  body  forces 

{p}   the  external  surfaces  pressures 

{F  }   any  external  force  applied  at  the  node 

{u}   displacements  at  a  typical  point 

{u  }   displacements  at  the  node  points 
n      r  r 

In  the  finite  element  method,  the  continuum  is  divided 
by  imaginary  lines  or  boundaries  into  a  finite  number  of 
small  elements.   These  "finite-elements"  are  assumed  to  be 
interconnected  at  a  discrete  number  of  boundary  points, 
called  nodes.   The  displacements  of  the  node  points  are 
interpreted  as  the  actual  displacements  of  the  corresponding 
points  of  the  soil  and  are  the  basic  unknowns  of  the  problem. 
In  the  single  phase  procedure  assumptions  are  made  relating 
the  displacements  within  an  element  to  the  nodal  displacements, 
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then  with  the  soil  constitutive  law,  substituted  into  the 
equilibrium  equation  reducing  it  to  the  matrix  form 

{P}e  =  [K]  {un}e  (5-8) 

where 

{P}e   is  the  vector  of  known  nodal  forces 

{u  }    is  the  vector  of  unknown  nodal  displacements 
n  r 

[K]   is  the  stiffness  matrix  relating  the  above. 
The  superscript  (e)  indicates  one  particular  element. 

This  is  not  however  possible  in  the  two-phase  problem 
since  in  the  equilibrium  equation  (5-7)  the  first  term  on 
the  right  hand  side  contains  {S}  the  unknown  pore  pressure 
which  cannot  be  expressed  in  terms  of  {u  }.   Equation  (5-7) 
can  be  reduced  to  matrix  form  in  the  same  manner  as  for  the 
single  phase  but  the  resulting  equation  has  an  extra  term. 

{P}e  -  [K]  {un}e  +  [K']  (H)e  (5-9) 


where 


(H)e  =  S/Kb 

S   is  the  unknown  pore  pressure  assumed  constant 
throughout  the  element 
K,    is  the  soil  bulk  modulus 

Another  equation  is  now  required  in  order  to  solve  for 
the  additional  unknown  resulting  from  the  introduction  of 

the  second  phase.   This  equation  equates  the  volumetric 

AY 
strain  (e.)  and  fluid  outflow  -£-   during  consolidation. 


from  which  we  obtain 


AV   =    fe    1 


{T}e    =    [K']T    {un)e 


where 


and 


{T}' 


,iT 


Xj    Kb  Area 


[K']x    {un}    -    (e^)    Kb   Area 
Equations    (5-9)    and    (5-11)    can  be    combined   to    form 
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(5-10) 


(5-11) 


(5-12) 


(5-13) 


K 


K 


T 
K'1    0 


(5-14) 


where 


and 


is  the  vector  of  known  loads  and  volumetric  strains 
fu 


n 


H 


is  the  vector  of  unknown  nodal  displacements 


and  pore  pressures.   The  equation  formulation  is  demonstrated 
below.   The  displacement  function  chosen  for  constant  strain 
triangular  elements  is 


u  =  ctj  +  a2  x  +  a3  y 
v  =  ak   +  a5  x  +  a6  y 


(5-15) 


or 


{u}  =  j    |  =  [M]  {a} 


(5-16) 
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The  strains  are  defined  as 


{e}  = 


CX 

y 

Y 
'xy 

6u 
37 

6v 
5y 

■  = 

sy     ^x 

0  10  0  0  0 
0  0  0  0  0  1 
0   0   10   10 


Oil 
0(2 

as 
a  if 
a  5 

a6 


(5-17) 


or 


{e}  =  [L]  {a} 


The  volumetric  strain  is 


(5-18) 


(e  )  =  (ex  +  ey)  = 


6u   6v 
6x   6y 


=  [010001]  {a}  ■  [g]  {a}  (5-19) 


Substituting  the  known  node  point  coordinates  into  equa 
tion  (5-15)  gives  the  relationship  between  node  point  dis- 
placements and  the  {a}'s. 


u. 

1 

vi 

uj 

.  = 

v. 
3 

uk 

.  V 

1  x.   y . 
i  'i 

0  0   0 

ooo 

1  xk  ^k 

0   0    0 


xi  y± 


xj  yj 


0   0 

xk  ^k 


(  Ol 

a2 

a3 

■ 

««, 

a5 

a6 

(5-20) 


or 


{un}  «  [A]  {a} 


(5-21) 
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Inverting  equation  (5-21) 

{a}  =  [A]"1  {u  } 


(5-22) 


Substituting  equation  (5-22)  into  (5-16),  (5-18)  and  (5-19) 
gives 


{u}  =  [M]  [A]"1  {un}  =  [N]  {un> 


{e}  =  [L]  [A]'1  {un>  =  [B]  {un} 


(O  =  [g]  [A]"1  {un}  =  [G]  {un} 


(5-23) 
(5-24) 
(5-25) 


The  elastic  stress-strain  relationship  is 


{a}  =  [D]  {e} 


(5-26) 


Substituting  equations  (5-23),  (5-24),  (5-25)  and  (5-26) 
into  the  two  equilibrium  equations  (5-7)  and  (5-10)  we  obtain 


[B]T  [D]  [B]  dv  {un}e  + 


[G]T  Kb  dv  (H) 


■I. 


[B]1  [D]  [eo]  dv  * 


f  {b}T  [N]  dv 


{p}T  [N]  ds  +  {F}^ 


(5-27) 


and 


^  ■  [G]  (un)e 


C5-28) 
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These  can  be  combined  as  before  to  give 


1  e 
P 


r         K 

0 


K'T 


(5-29) 


where 


{P}6   =    {F}*   +    f    [B]T    [D]    {eo}    dv   + 


+    f    {p}T    [NJ    ds 
'  s 

{T}e  =  [^)  (K"b  Area) 
[K]  =  f  [B]T  [D]  [B]  dv 


{b}1  [N]  dv 


=  1 


[K»]  =    [G]1  dv  K, 


[K'T]  =  f  [G]  dv  K, 
'  v 


(5-30) 

(5-31) 
(5-32) 
(5-33) 
(5-34) 


{u  }   =  the  unknown  modal  displacement 
(H)   =  the  unknown  pore  pressure. 

[B] ,  [D] ,  [G] ,  [N]  and  thus  [K]  and  [K']  are  functions  of 

location. 

(b>,  (p),  {F  },  {eQ}  anc*  {T}  are  loads  which  are  to  be 

specified. 

{u  }   and  (H)   can  hence  be  solved. 
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5.3   The  Consolidation  Process 


The  solution  to  the  transient  flow,  or  consolidation 
problem,  follows  directly  from  the  two  phase  formulation  by 
specifying  values  for  the  volume  changes  in  the  load  vector. 


Assume 


1  AV 
At  T 


6v: 

6x 


<5v 


6y  j 


and 


Darcy's  law 


vx  = 


"xx  2Tx 


(5-35) 


(5-36) 


v   -  -  k    6h 

vy    kyy  *y 


Substituting  equations  (5-36)  ,  and  h  =  S/y   into  equa- 


tion   (5-35)    gives 


1_ 

At 


'AVl         Jjcx   62S        ^XZ  62s 
I  V  J    "   Yw      GF        yw      6y-T 


or 


,       K,  •  Area  r- 
{T}  =  -A- 


AT 


62S  +     62S 
xx  6x2    yy  6yz 


(5-37) 


(5-38) 


The  pore  pressure  S  has  been  defined  as  constant  throughout 
each  element.   If  quadrilateral  elements  are  used  in  the  finite 
element  program  then  a  quadrilateral  and  its  four  adjacent 
quadrilaterals  allow  the  $'s  in  the  following  polynomial 
to  be  defined 


S  =  B,  +  B2  x  +  33 


y  +  B  x2  +  $,  y: 


(5-39) 
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W" 


=    [E] 


(5-40) 


.*.       {$}    =    [E]"1    {Sn> 


From  equation    (5-39) 
62S 


JyTT    _     235 


53^  "  2e<.  ' 
From  the  relevant  portion  of  equation  (5-41) 


E|»l         Ei,2         Ei»3         Sit         Ei*5 
E5  1         E5  2         E5  3         E5k         E5  5 


S. 

1 


Substituting   equations    (5-42)    and    (5-43)    into    (5-38) 
Area    r 

—  C 


lC,    Area 
{T}    -    At    •    \—     |, 


(5-41) 


(5-42) 


(5-43) 


Ekl    S.    +   E\2    S..    +    £  L,.,.      • 


f  3    "k         ^"m       ^ 


+     Fuc     S. 


+  k. 


yy 


E5l    S.    ♦   E52    Sj    *   E53    Sk   + 


EKU    S0    +   Ecc    S 

5«t         £  5  5m 


(5-44) 
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where 

{T}  is  the  volume  change  for  each  element 

At  is  the  specified  time  increment 

K,  is  the  soil  Bulk  modulus 

Area  is  the  element  area 

k   ,  k    the  principle  permeabilities. 
xx'   yy      y  v       r 

The  solution  procedure  for  the  consolidation  problem 
consists  of  alternately  solving  equations  (5-14)  and  (5-44) 


K 

K'  " 

T 
K*1 

0 

■ 

(5-14) 


K,  Area 

{T>  =  At  -2- [ff]  {S  } 

Tw 


(5-44) 


First,  the  initial  undrained  displacements  and  pore  pressures 
are  found  by  assuming  {T} ,  the  volume  changes,  to  be  zero  in 
equation  (5-14)  .   These  pore  pressures  are  then  substituted 
into  equation  (5-44)  and  new  volume  changes  calculated  after 
dissipation  for  some  time  increment  At.   The  {T}'s  are  now 
used  as  known  volume  changes  in  the  load  vector  in  (5-14) 
to  determine  new  displacements  and  pore  pressures.   This 
cycle  is  repeated  until  the  pore  pressures  have  dissipated 
to  some  specified  level. 
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5.4   The  Finite-element  Program 
The  original  finite  element  program  CONSOL  consisted  of 
a  MAIN  program  and  ten  major  subroutines.   Four  supplementary 
subroutines  can  be  employed  for  plotting,  special  consolida- 
tion output,  rigid  loads  and  time  dependent  loading.   A 
further  subroutine  QUASI  has  been  added  to  take  into  account 
the  special  soil  constitutive  behavior  and  will  be  described 
later  in  detail.   The  major  program  divisions  are  listed 
below  with  a  brief  description  of  their  function.   Input 
formats  and  a  program  listing  are  given  in  Appendix  D. 

MAIN  -  routes  the  finite  element  procedure  through  the  neces- 
sary subroutines  for  a  solution  to  the  input  problem. 

MLDATA  -  reads  in  material  properties  and  calculates  all 
necessary  soil  parameters . 

DATA1  -  reads  in  node  point,  element,  pore  pressure  and 
boundary  loading  data.   Sets  up  location  arrays  and 
determines  the  band  width. 

DATA2  -  reads  in  boundary  drainage  conditions,  determines 

element  incidences  and  terms  for  the  flow  calculations. 

SOLSTF  -  formulates  the  load-volumetric  strain  vector  and 
the  overall  stiffness  matrix. 

MODIFY  -  takes  into  account  specified  displacements  and  pore 
pressures . 

QUAD  -  assembles  the  stiffness  matrix  for  the  quadrilateral 
element  by  subdividing  it  into  four  triangles  and  then 
superimposing  results. 


81 

BANSOL  -  triangularizes  the  stiffness  matrix  by  Gaussian 

Elimination  and  solves  for  the  unknown  displacements 

and  pore  pressures  by  backsubstitution. 
FLOWCA  -  calculates  the  volumetric  strains  due  to  each 

increment  of  fluid  flow  during  consolidation. 
TMCRIT  -  calculates  the  time  increments. 
STRESS  -  from  the  determined  displacements  calculates 

strains  and  stresses  and  prints  out  results  in  tabular 

form. 
QUASI  -  checks  if  the  stress  path  has  crossed  the  quasi- 

surface  and  if  so  changes  the  element  soil  properties. 

Two  major  changes  were  made  to  the  MIT  CONSOL  program 
and  are  described  below.   In  the  original  program  input  data 
were  read  in  the  three  DATA  subroutines,  the  stiffness 
matrix  formulated  in  SOLSTF  and  the  initial  undrained  results 
calculated  in  BANSOL  and  output  in  STRESS.   For  the  consoli- 
dation phase  the  program  then  entered  a  loop  consisting  of 
BANSOL  (back-substitution  only),  FLOWCA  and,  whenever  output 
was  desired,  STRESS.   Each  cycle  through  this  loop  represents 
one  time  increment  of  excess  pore  pressure  dissipation. 

A  major  revision  of  the  existing  program  was  necessary 
to  satisfy  the  requirements  of  the  present  research.   For  the 
soil  model  being  considered  a  check  of  the  effective  stress 
state  against  the  q-yield  surface  is  required  after  each 
dissipation  increment.   The  program  was  therefore  altered 
such  that  subroutine  STRESS  was  called  during  every  cycle 


82 

but  that  the  write-out  section  within  STRESS  was  again  only 
utilized  at  specified  times. 

The  second  major  change  made  to  the  original  program 
was  the  addition  of  the  soil  model  subroutine,  QUASI.   This 
is  called  from  STRESS,  for  each  grid  element,  after  the  dis- 
placements, strains  and  stresses  have  been  calculated  at  a 
particular  time.   The  insitu  effective  stresses  are  deter- 
mined in  QUASI  assuming  a  simple  geostatic  state  with  water 
table  at  the  soil  surface  and  with  assumed  values  of  soil 
unit  weight  and  coefficient  of  earth  pressure  at  rest. 
These  insitu  stresses  prior  to  loading  are  required  both 
in  the  calculation  of  the  octahedral  shear  stress  increment, 
At  ct,  and  as  a  representation  of  the  field  q-consolidation 
pressure.   For  this  pressure  a  corresponding  q-yield  curve 
is  determined  and  a  comparison  made  with  the  element  effective 
octahedral  stress  increments  due  to  the  loading.   If  the 
incremental  stress  state  represents  a  point  outside  the 
elliptic  yield  curve  then  the  post  yield  elastic  parameters 
are  assigned  to  that  soil  element  and  are  used  in  the  next 
cycle  determination.   This  change  in  element  soil  parameters 
necessitates  leaving  the  dissipation  loop  described  above 
and  calling  SOLSTF  in  order  to  calculate  a  new  stiffness 
matrix.   Triangularization  is  then  required  within  BANSOL. 
These  re~routings  have  been  inserted  into  the  program  ^and  a 
simplified  flow  chart  is  shown  in  Figure  5.1. 
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5.  5   Illustrative  Problems 

5.5.1   Introduction.   In  order  to  determine  the  effect 
of  the  quasi-preconsolidation  model  on  the  consolidation 
process,  finite  clement  solutions  were  obtained  for  four 
different  consolidation  problems.   These  represented  the 
one-  and  two-dimensional  cases  shown  in  Figure  5.2,  with 
and  without  consideration  of  the  quasi-preconsolidation 
behavior.   Plots  and  discussion  of  these  results  are  given 
in  subsequent  sections. 

In  the  finite  element  method  the  degree  of  freedom  of 
each  node  point  must  be  specified.   A  choice  exists  as  to 
this  freedom  for  the  bottom  row  of  node  points,  representing 
the  soil-rock  interface.   Each  node  may  be  fixed,  or  may  be 
permitted  to  move  only  in  the  x-direction.   These  cases 
represent,  respectively,  the  u  =  0  (rough]  and  x   =  0  (smooth) 
boundary  conditions  of  analytical  stress  distribution  theory. 
Finite  element  programs  for  the  two  cases  were  run,  and 
initial  excess  pore  pressures  calculated  using  the  equation 

Au  =  Ao\  .  +  (0.2)  At   .  (5-45) 

oct  oct 

Contours  of  excess  initial  pore  pressure  are  plotted  in 
Figures  5.3  and  5.4.   When  the  lower  nodes  are  allowed  a 
single  degree  of  freedom,  the  fixed  right  hand  grid  boundary 
is  the  only  component  preventing  lateral  deformation. 
Horizontal  stresses  do  not  dissipate  to  zero  with  increasing 
distance  from  the  load,  as  reflected  in  the  pore  pressure 
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Uniformly  Distributed   Load 


i    i    i    i    i    i    i    i    i    i    i    r 


T 
Drainage 


Clay 


No  Drainage 

Rock 


a  )- One- Dimensional  Consolidation  Configuration 


Drainage 


Uniformly  Distributed   Strip   Load 

i    i'   i    i 


Clay 


No  Drainage 

Rock 
b)-  Two-  Dimensional    Consolidation   Configuration  . 


FIGURE    5.2     ONE     AND   TWO    DIMENSIONAL     CONSOLIDATION 
CONFIGURATIONS. 
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contours  of  Figure  5.4.  Thus,  it  appears  that  fixed  lower 
boundary  node  points  are  more  realistic.  It  was  therefore 
decided  to  fix  the  bottom  node  points  in  all  cases. 

5.5.2.   One -Dimensional  Consolidation  Problem,  No 
Quasi -Preconsolidat ion  Effect.   A  finite  element  grid, 
consisting  of  seventy-five  elements  and  shown  in  Figure  5.5, 
was  employed  for  the  one-dimensional  no  yielding  solution. 
In  order  to  obtain  one-dimensional  deformation  and  drainage 
behavior,  a  uniformly  distributed  load  was  assigned  over 
the  upper  surface,  with  drainage  permitted  only  from  the 
same  boundary.   This  problem  allows  a  comparison  between 
the  finite  element  solution  and  the  analytical  Terzaghi  one- 
dimensional  consolidation  theory.   For  the  case  where  the 
initial  excess  pore  pressure  u  =  u  is  constant,  the 
Terzaghi  solution  may  be  written 

u  -  uQ  I     1  sin  M   exp  -M2T  (5.46) 


where 


m  _  C  2m  +1 1  it  n   ,   ,   - 

M  =  ■» -z — *—       m=  0,  1 ,  L  ,  3  .... 


T 


c  t 
v 


and  H  -   thickness  of  the  clay  layer. 

A  simple  computer  program  was  developed  to  solve  this  equa- 
tion, calculating  u/uQ  for  input  values  of  depth  factor  Z/H 
and  time  factor  T.   The  program,  in  which  a  summation  of  m 
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from  0  to  50  was  employed,  is  listed  in  Figure  5.6.   In 
Figure  5.7  are  plotted  isochrones  for  different  values  of 
the  time  factor.   The  Terzaghi  solution  is  given  in  the 
upper  portion  of  the  figure,  and  that  for  the  five  row 
finite  element  program  in  the  lower  part.   Agreement  even 
for  this  coarse  grid  is  quite  good,  except  near  the  drainage 
boundary  where  pore  pressures  in  the  finite  element  method 
dissipate  more  rapidly  than  in  the  Terzaghi  solution. 
.Christian  and  Boehmer  (1970) ,  compared  solutions  from  the 
original  CONSOL  program  for  various  numbers  and  shapes  of 
grid  elements,  with  Terzaghi's  solution,  finite  difference 
solutions,  and  an  analog  computer  solution.   For  a  ten  row 
finite  element  solution  plots  were  indistinguishable  from 
Terzaghi's  solution.   For  a  five  row  program  they  also  found 
some  deviation  from  the  correct  solution  but  that  the  iso- 
chrones throughout  most  of  the  depth  were  close  to  the 
correct  values. 

In  Figures  5.8,  5.9  and  5.10  are  plotted  excess  pore 
pressure  ratio  values  versus  logarithm  of  the  time  factor 
for  elements  in  the  top,  second  and  third  rows  of  the  finite 
element  grid.   The  Terzaghi  solutions  are  also  plotted  and 
illustrate  the  close  agreement  except  in  the  first  row  of 
elements  adjacent  to  the  drainage  boundary.   Surface  settle- 
ments are  plotted  in  dimensionless  form  in  Figure  5.11. 
These  relative  settlements  are  the  actual  settlements  p 
divided  by  the  ultimate  consolidation  settlement  Puit« 
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DIMENSION  Z(6) ,T(10) 

PI=3. 14159265 

READ(5,1)   (T(K) ,  K=l,10) 

1  FORMAT (8F10.0) 
DO  4    J=l,5 
AJ  =  J 

Z(J)  =  (2.0*AJ  -  1.0)  /  10.0 

DO  4    K=l,10 

SUM  =0.0 

DO  2    1=1,51 

M  =  I  -  1 

AM  =  M 

BIGM  =  (2.0*AM  +  1.0)  *  PI  /  2.0 

Fl  =  BIGM  *  Z(J) 

S  =  SIN(Fl) 

F2  =-BIGM  *  BIGM  *  T(K) 

E  =  EXP(F2) 

F3  =  (2.0  /  BIGM)  *  S  *  E 

2  SUM  =  SUM  +  F3 
WRITE(6,3)   Z(J) ,  T(K) ,  SUM 

3  FORMAT(  /7HDEPTH  = ,F8 . 5 , 5X, 3HT  = , F9 . 5 , 5X, 3HU  =,F6.3/  ) 

4  CONTINUE 
END 

FIGURE    5.6     ANALYTICAL     CONSOLIDATION      PROGRAM. 
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FIGURE    5.7-  PORE     PRESSURE     ISOCHRONES     FOR 
ANALYTICAL    AND    FINITE    ELEMENT 
SOLUTIONS. 
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In  Section  5.5.3  these  plots  will  be  compared  with  results 
obtained  for  a  similar  grid  configuration,  but  with  con- 
sideration of  the  quasi-preconsolidation  behavior. 

5.5.3   One -Dimensional  Consolidation  Problem  with  Quasi- 
Preconsolidation  Effect.   A  finite  element  solution  was 
obtained  to  the  corresponding  one-dimensional  consolidation 
problem  in  which  the  quasi-preconsolidation  effects  previ- 
ously discussed  were  included.   This  solution  was  generated 
using  the  seventy-five  element  mesh  shown  in  Figure  5.5. 
In  addition,  in  order  to  observe  in  detail  changes  within 
an  element  when  it,  or  a  neighboring  element  yields,  two 
further  meshes  consisting  of  one  and  two  elements  were  used. 
These  are  discussed  below. 

The  single  element  representation  illustrates  the  effect 
which  sudden  change  in  elastic  modulus,  corresponding  to 
yielding,  has  on  excess  pore  pressures  and  volume  strain  at 
the  point  where  yield  occurs.   The  single  element  was  a 
square,  undrained  on  three  sides  and  subjected  to  a  uniform 
vertical  surface  pressure.   As  shown  in  Figure  5.12  the 
nodes  were  constrained  to  permit  only  one-dimensional  dis- 
placements of  the  element  boundaries.   Figures  5.13  and  5.14 
show,  respectively,  plots  of  pore  pressure  ratio  and  relative 
settlement  versus  logarithm  of  the  time  factor. 

The  effect  of  yielding  was  obtained  by  reducing  the 
elastic  modulus  E  by  a  factor  of  four  after  five  time 
increments  of  excess  pore  pressure  dissipation.   This 
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FIGURE     5.12      SINGLE      ELEMENT    CONSOLIDATION     GRID. 
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induced  yielding  resulted  in  the  pore  pressure  and  relative 
settlement  response  depicted  respectively  in  Figures  5.15 
and  5.16.   The  curves  for  the  non  yielding  case  are  plotted 
for  comparison.   The  sudden  yielding  and  increase  in  com- 
pressibility of  the  soil  skeleton  results  in  a  transfer  of 
load  to  the  pore  fluid,  and  the  consequent  increase  in  excess 
pore  pressure  ratio  from  0.76  to  0.94  (Figure  5.15).   For  the 
single  element  representation,  these  values  can  be  verified 
directly:   For  the  one-dimensional  consolidation  process, 
applied  boundary  loads  are  supported  initially  by  the  pore 
water  and  not  by  the  soil  skeleton.   The  initial  stresses, 
normalized  with  respect  to  the  applied  load,  are 

Total  stress  a   =  1.00 

Excess  pore  pressure     u  =  1.00        (5-47) 

Effective  stress         a   =  0.00 

With  time,  as  water  drains  from  the  clay  soil  the  excess 
pore  pressure  at  a  point  decreases  and  the  effective  stress 
increases.   The  one  element  representation  indicates  a  pore 
pressure  ratio  value  of  0.76  just  prior  to  yielding.   Since 
in  the  one-dimensional  case  the  total  stress  remains 
constant  the  normalized  stresses  at  this  time  are 


Total  stress  a  =  1.00 

b 

Excess  pore  pressure    u^  =  0.76        (5-48) 

Effective  stress        a,  =  0.24 
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where  the  subscripts  b  and  a,  respectively  denote  behavior 
just  before  and  just  after  yield. 

For  a  model  assuming  elastic  skeleton  behavior  the 
effective  stress  and  the  strain  are  related  by  Hookes  Law 


E  =  £ 


(5-49) 


and  the  effective  stress  principle  for  the  state  just  prior 
to  yielding  may  be  written 


from  which 


°b  =  Ub  +  Eb  '  £ 


e  =  J> b 


(5-50) 


(5-51) 


Just  after  the  instantaneous  yielding  the  strain  has  the 
same  value  as  just  before,  since  no  flow  of  water  out  of 
the  soil  has  taken  place. 


Therefore 


ua  =  aa  -  Ea 


ab  •  ub 


:b   J 


Substituting  in  the  known  value 


u   =  1.00  -  0.25  (1.00  -  0.76) 


=  0.94 


(5-52) 


=  the  pore  pressure  calculated  in  the  program 
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In  Figure  5.16  relative  settlements  are  plotted  with 
time  for  the  no-yielding  and  with  yielding  conditions.   The 
decrease  in  elastic  modulus  upon  yielding  reduces  the  coef- 
ficient of  consolidation  c  ,  and  consequently  the  rate  of 
pore  pressure  dissipation.   Larger  settlements  also  result 
from  the  reduced  E. 

The  results  for  the  two  element  grid  shown  in  Figure 
5.17,  illustrate  the  effects  which  yielding  in  one  element 
has  on  a  neighboring  element.   Using  this  grid  configuration, 
two  features  can  be  observed  clearly:   the  effect  of  yielding 
of  the  upper  element  on  the  lower  element  and,  that  of  the 
subsequent  yielding  of  the  lower  element  on  the  already 
yielded  upper  element.   When  no  yield  occurs,  pore  pressure 
dissipation  and  relative  settlement  for  each  element  are  as 
shown  in  Figures  5.18  and  5.19  respectively.   The  pore 
pressure  in  Element  1  initially  dissipates  more  rapidly 
than  in  Element  2  because  the  only  drainage  surface  is  the 
top  boundary.   At  later  times  however,  its  drainage  is 
delayed  as  water  from  the  lower  element  drains  through  it. 
Thus,  the  pore  pressure  approaches  zero  in  both  elements  at 
the  same  time.   In  Figure  5.19  are  plotted  relative  settle- 
ments of  the  elements.   These  are  obtained  by  dividing  the 
vertical  displacements  of  Node  Point  1  which  represent  the 
total  surface  settlement,  the  vertical  displacements  of 
Node  Point  2  which  represent  the  settlement  within  Element  2, 
and  the  difference  which  represents  the  settlement  within 
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Element  1,  by  the  ultimate  settlement  of  Node  Point  1.   As 
in  the  case  of  the  pore  pressure  dissipation,  settlement 
occurs  initially  more  rapidly  in  the  upper  element  adjacent 
to  the  drainage  boundary.   Since  the  elements  are  identical 
final  settlements  within  each  are  equal  and  one  half  of  the 
total  surface  settlement. 

The  effect  of  yielding  was  observed  by  reducing  instan- 
taneously the  elastic  modulus  by  a  factor  of  four  in  the 
upper  element  after  five  dissipation  cycles  (at  T  =  0.0075) 
and  in  the  lower  element  after  fifteen  cycles  (at  T  =  0.0186) 
The  dimensionless  time  factor  T  is  defined  by 

cv 
T  =  jjr  •  t  (5-53) 

in  which  c  refers  to  the  now  yielded  coefficient  of  con- 
solidation.  For  an  elastic  analysis 

cy  =  £  •  D  (5-54) 

'w 

where  D  -  2G  ^X=  (5-55) 

l-2v 

G  =  |  •  (l+2v)  (5-56) 


and 


2 


c  =  coefficient  of  consolidation 
v 

t  =  real  time 

H  =  maximum  drainage  path  length 

k  =  coefficient  of  permeability 

Yw  =  unit  weight  of  water 
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D  =  constraint  modulus 

G  =  shear  modulus 

v  =  Poisson's  ratio  of  the  soil  skeleton 

E  =  Young's  modulus 

The  excess  pore  pressure  ratio  versus  logarithm  of  the 
time  factor  for  Element  1  is  given  in  Figure  5.20,  along  with 
the  no  yield  curve.   When  yielding  occurred  in  Element  1 
there  was  again,  as  in  the  one  element  program,  an  instan- 
taneous increase  in  excess  pore  pressure  due  to  the  increased 
compressibility  of  the  skeleton  and  the  transfer  of  the  load 
to  the  pore  water.   When  Element  2  yielded  at  time  factor 
0.0186  ,  the  effect  of  its  pore  pressure  increase  on  Element 
1  was  small.   The  pressure  gradient  due  to  the  drainage 
boundary  dominated  behavior  in  Element  1,  and  the  effect  of 
the  increased  gradient  and  increased  flow  of  water  from 
Element  2  upon  its  yielding  was  insignificant.   The  effect 
of  the  yielding  of  an  element  on  another  element  which  has 
already  yielded,  but  which  does  not  bound  a  drainage  zone, 
will  be  observed  in  the  results  for  the  seventy-five  element 
mesh  discussed  below. 

Pore  pressure  behavior  for  the  lower  element,  Number  2, 
is  plotted  in  Figure  5.21,  along  with  its  no  yield  curve. 
When  the  upper  element  yielded  at  T  =  0.0075  its  pore 
pressure  increased  instantaneously,  reducing  the  gradient 
and  flow  out  of  Element  2.   The  rate  of  dissipation  is 
markedly  decreased,  and  the  plot  is  displaced  to  the  right 
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of  the  no  yield  curve.   At  time  factor  0.0186  when  it 
yielded,  the  pore  pressure  in  Element  2  increased  from  88 
to  97  percent  of  its  initial  value,  then  dissipated  w.ith 
time  to  zero. 

Relative  settlements  within  the  upper  and  lower  elements 
and  the  total  relative  settlement  for  the  yielding  case, 
along  with  the  corresponding  no  yield  curves,  are  given  in 
Figures  5.22,  5.23  and  5.24,  respectively.   Settlement  due 
to  Element  1,  Figure  5.22,  is  similar  in  form  to  that  found 
for  the  single  element  mesh.   After  yielding,  the  rate  of 
settlement  is  decreased  and  the  final  magnitude  increased. 
No  effect  is  observed  on  this  plot  when  Element  2  yields. 
In  the  lower  element,  Figure  5.23,  settlements  initially  are 
equal  to  the  no  yield  values.   After  Element  1  yields,  dis- 
placements are  small  until  Element  2  itself  yields  at  time 
factor  0.0186.   Settlements  then  occur  in  greater  magnitude 
than  for  the  no  yield  behavior.   In  Figure  5.24  is  given  the 
total  relative  surface  settlement  representing  the  sum  of 
the  previous  two  plots,  and  showing  the  increase  in  settlement 
due  to  the  increased  compressibility  of  the  soil  skeleton 
upon  yielding. 

The  seventy- five  element  mesh  was  used  with  the  labora- 
tory derived  model  incorporated  into  the  program  to  determine 
solutions  to  the  one-dimensional  consolidation  problem  for  a 
soil  which  exhibits  a  quasi-preconsolidation  effect.   Rather 
than  inducing  yielding  at  a  specific  time,  as  in  the 
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illustrative  one  and  two  element  representations,  yielding 
occurred  when  the  effective  stress  state  in  an  element 
exceeded  the  elliptic  cone  yield  surface  defined  in  Section 
4.3.   After  each  cycle  of  pore  pressure  dissipation  the 
effective  octahedral  stress  increments  Aa   .  and  La or*»  for 
each  element,  were  checked  to  determine  whether  they  lay 
within  or  outside  the  effective  stress  yield  surface  defined 
by  the  equation 


(AToct  >   =  0.032818  (o  )z  +  0.067384  (a  )  (Aooct  ) 
q  M  H  q 

-  0.52570  (Aa    )2  (4-7) 

octqJ 


where 


At   .   is  the  octahedral  shear  stress  increment  at  yield 

q 

Act  „«.   is  the  octahedral  effective  normal  stress  incre- 
oct 

q 

ment  at  yield 
a      is  the  q-consolidation  pressure 

This  was  accomplished  in  the  following  steps: 

1.  For  the  particular  element  the  effective  octahedral 
normal  stress  at  its  center  was  calculated.   This 
represented  the  q-conf inement  pressure  for  that 
element . 

2.  For  this  q-conf inement  pressure  the  yield  curve 
on  the  At    -  Aa t  plane  was  determined 
(Figure  5. 25) . 
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FIGURE  5.25      ELEMENT     YIELD     AND     NO    YIELD    CONDITIONS. 
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3.  For  the  element  Aa   .  the  Ax   .  was  compared  with 

oct       oct 

the  ordinate  of  the  curve  intersection  Ax 

oct 

q 

4.  If  Ax    was  greater  than  At   t  the  element  had 

q 
yielded,  and  soil  properties  were  reassigned.   If 

AxQct  was  less  than  Ax     the  element  had  not 

q 

yielded,  and  soil  properties  remained  at  pre-yield 

values. 
Sketches  of  no  yield  and  yield  conditions  are  shown  in 
Figure  5.25. 

The  state  of  the  soil  relative  to  yielding  at  each 
depth  for  the  seventy-five  element  representation  is  shown 
schematically  in  Figure  5.26.   At  time  zero  when  the  surface 
load  is  applied,  the  excess  pore  pressure  ratios  are  equal 
to  unity  in  all  elements  for  the  one-dimensional  case,  and 
none  of  the  elements  have  yielded.   This  is  illustrated  in 
the  left  hand  diagram  in  Figure  5.26  by  the  unshaded  squares. 
Due  to  the  hydraulic  gradient  across  the  drainage  boundary, 
water  flow  out  of  the  soil,  and  the  excess  pore  pressures 
dissipate  with  time.   Since  both  loading  and  drainage  are 
one-dimensional,  behavior  is  identical,  within  the  accuracy 
of  the  numerical  method,  in  all  elements  in  any  row  of  the 
finite  element  grid.   Pore  water  drainage  occurs  most 
rapidly  in  the  top  row  of  elements  adjacent  to  the  drainage 
boundary,  and  at  time  factor  0.00046  the  excess  pore  pressure 
ratio  has  dissipated  to  0.9569.   The  effective  stresses  have 
changed  such  that  at  this  time  the  effective  octahedral 


]  2  0 


bo 

si  ii 


CD 


^^^ 

ro 


ill 

1 

do 

ii  H 


do 

H-3 


O 

»oo 

OC0 

do 

u      II 
J-   3 


CD 

«frCT> 

OCD 

om 

OCT) 

do 

h    n 
h-  3 


o 

2 


O^j 

0     u 

h"   3 


o 

5 

q: 

id 

«r 

=> 

oo 
en 

(J 

cc 

Q. 

ID 
K 
O 
Q. 

oo 

CO 
ID 

X 

ID 

D 

-J 
UJ 

>-    • 

O 

<H 
2 

a:  id 
oS 
•-H 

gUj 

i=£ 

(M 

tf> 

ID 

a: 

D 
O 

LI 


121 

stress  state  lies  outside  the  yield  surface,  Figure  5.26. 
This  result  is  indicated  schematically  in  the  second  diagram 
in  Figure  5.26.   Thus  the  elastic  modulus  was  reduced  in 
all  top  row  elements  to  the  post  yield  value.   At  T  =  0.03050, 
when  the  excess  pore  pressure  ratio  in  the  second  row  of 
elements  has  dissipated  to  0.8780,  the  second  row  yields. 
The  block  diagrams  in  Figure  5.26  illustrate  the  time 
factor  at.  which  this  and  lower  rows  yield,  and  the  excess 
pore  pressure  ratios  within  the  yielding  elements. 

In  Figure  5.27  is  shown  the  excess  pore  pressure  ratio 
versus  logarithm  of  the  time  factor  for  a  top  row  element 
for  the  seventy-five  element  one-dimensional  consolidation 
representation,  with  the  quasi-preconsolidation  model 
incorporated.   Also  plotted  for  comparison  is  the  no  yield 
curve  from  Figure  5.8.   As  in  the  simple  one  and  two  element 
cases,  excess  pore  pressures  increase  suddenly  when  the 
element  yields.   And,  as  in  the  two  element  model,  there  is 
no  significant  effect  on  these  pore  pressures  when  lower 
elements  subsequently  yield. 

A  similar  plot  for  a  second  row  element  in  Figure  5.28 
shows  four  distinct  effects  due  to  yielding  of  the  top, 
second,  third  and  fourth  rows.   The  top  row  of  elements 
yielded  at  time  factor  0.00046,  off  the  scale  of  this  plot. 
The  effect  of  this  was,  as  described  before,  to  reduce  both 
the  gradient  and  flow  into  this  upper  region.   The  excess 
pore  pressure  ratio  plot  is  therefore  displaced  to  the  right 
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of  the  no  yield  curve.   At  time  factor  0.03050  row  two 
yields  and  the  load  transfer  caused  the  excess  pore  pressure 
increase  previously  described.   Dissipation  of  pore  water 
pressure  then  progresses  smoothly  until  T  =  0.16982,  when 
the  third  row  yields.   In  this  case,  unlike  the  top  row  of 
elements,  the  excess  pore  pressures  are  effected  by  the 
yielding  of  elements  directly  below.   The  effect  of  the 
increased  flow  of  water  into  the  element  is  significant 
resulting  in  a  further  slowing  of  the  consolidation  process. 
This  is  observed  on  the  plot  as  another  displacement  of  the 
excess  pore  pressure  curve  to  the  right.   A  smaller  effect 
can  be  noted  at  time  factor  0.44384  when  the  fourth  row 
yields. 

In  Figure  5.29  the  plots  for  a  third  row  element  are 
given.   The  top  row  has  already  yielded  by  time  factor  0.10, 
and  the  curve  is  offset  from  the  no  yield  behavior.   At  time 
factor  0.03050  row  two  yields,  resulting  in  a  curvature 
change  at  this  time  and  a  further  retarding  of  the  rate  of 
pore  pressure  dissipation.   At  T  =  0.16982  the  element  under 
consideration  yields  with  the  associated  sudden  pore  pressure 
increase.   The  fourth  row  yields  at  time  factor  0.44384  and 
the  fifth  row  yields  at  time  factor  0.94840.   In  each  case, 
the  effect  of  the  increased  compressibility  of  lower  rows  is 
to  reduce  the  rate  at  which  the  excess  pore  pressure  dissi- 
pates.  Similar  results  were  obtained  for  lower  rows  of 
elements . 
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Relative  surface  settlements  are  plotted  versus  log- 
arithm of  the  time  factor  in  Figure  5.30,  with  the  no  yield 
case  shown  for  comparison.   Yielding  of  the  top  row  of  ele- 
ments has  already  occurred  by  time  factor  .001  and  the 
settlement  curve  lies  to  the  left  of  that  for  no  yield 
behavior.   Changes  in  curvature  resulting  in  increased 
settlements  can  be  observed  as  each  subsequent  row  of 
elements  yields,  at  time  factors  0.03050,  0.16982  and 
0.94840.   in  order  to  better  illustrate  these  curvature 
changes  in  the  settlement  plot,  four  programs  were  run 
in  which  pore  pressures  corresponding  to  states  just  prior 
to  yielding  were  input.   Pore  pressure  dissipation  was  then 
continued  but  yielding  was  not  permitted.   These  results 
are  shown  in  Figure  5.30  as  dotted  curves,  and  illustrate 
the  deviation  in  the  settlement  curve  due  to  yielding 

5.5.4   Two -Dimensional  Consolidation  Problem,  No 
Quasi -Preconsolidation  Effect .   The  two-dimensional  problem 
was  introduced  in  Section  5.5.1  and  illustrated  in  Figure 
5.2.   It  consists  of  a  stratum  of  clay  soil  underlain  by  a 
rigid  impervious  boundary,  and  loaded  by  a  uniformly  dis- 
tributed surface  strip  load.   This  plane  strain  configura- 
tion is  represented  by  the  finite  element  grid  consisting 
of  100  square  elements  and  126  nodal  points,  shown  in 
Figure  5.31.   Output  from  the  two-dimensional  finite  ele- 
ment consolidation  program  consists  of  horizontal  and 
vertical  nodal  displacements  and  excess  pore  pressure 
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ratios  at  different  values  of  the  time  factor.   The  time 
factor  T  is  defined  as  in  Equation  (5-53)  with  II  equal  to 
the  thickness  of  the  clay  stratum.   The  excess  pore  pressure 
ratio  is  the  actual  pore  pressure  at  the  center  of  an  ele- 
ment, divided  by  the  magnitude  of  the  distributed  load. 

Contours  of  the  initial  distribution  of  excess  pore 
pressure  are  given  in  the  upper  diagram  of  Figure  5.32. 
Values  are  greatest  directly  under  the  center  of  the  load, 
and  decrease  with  depth  and  with  horizontal  distance  from 
the  loaded  area.   The  lower  diagram  of  Figure  5.32,  and  the 
upper  diagram  of  Figure  5.33,  show  the  rapid  dissipation 
of  pore  pressures  in  the  top  row  of  elements  adjacent  to 
the  only  drainage  boundary.   With  further  time  pore  pres- 
sures dissipate  in  all  elements,  and  the  contour  lines  move 
down  through  the  soil  layer.   By  time  factor  0.5967]  the 
pore  pressure  ratios  in  all  elements  are  less  than  0.1. 
Figures  5.32,  5.33  and  5.34  illustrate  with  contour  diagrams 
the  progression  of  pore  pressure  dissipation. 

Plotted  below  each  pore  pressure  contour  diagram  is  the 
corresponding  surface  settlement.   Initially,  at  time  factor 
equal  to  zero,  there  is  settlement  under  the  load  and  heave 
outside  the  loaded  area.   Had  the  finite  element  grid 
extended  further  laterally,  this  heave  would  have  decreased 
to  zero  at  some  distance  from  the  load.   The  area  enclosed 
by  the  settlement  portion  of  this  plot  is  equal  to  the 
heave  area,  indicating  that  the  initial  settlement  was  due 
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to  distortion  of  the  soil  with  no  volume  change.   Consolida- 
tion settlements  occur  with  time,  and  settlement  under  the 
center  of  the  loaded  area  increases,  while  the  peripheral 
heave  decreases. 

5.5.5   Two -Dimensional  Consolidation  Problem,  with 
Quasi-Preconsolidation  Effect.   The  two-dimensional  consoli- 
dation solution  for  the  case  in  which  the  quasi-preconsolida- 
tion  effect  is  present  was  obtained  using  the  grid  illustrated 
in  Figure  5.31  and  employing  the  laboratory  derived  soil  model 
as  described  for  the  seventy-five  element  one-dimensional 
program  in  Section  5.5.3.   Results  are  presented  in  the  form 
of  excess  pore  pressure  ratio  contour  diagrams  and  correspond- 
ing surface  settlements  at  different  time  factors. 

When  the  load  is  applied  at  time  zero  the  excess  pore 
pressure  distribution,  given  in  the  top  diagram  of  Figure 
5.35,  is  identical  to  the  no  yield  case  of  Figure  5.32. 
However  several  elements  in  the  two  rows  under  the  load 
yield  immediately  and  pore  pressures  increase  to  the  values 
shown  in  the  lower  diagram  of  Figure  5.35.   The  elastic 
modulus  for  a  yielding  element  is  reduced  by  a  factor  of 
four,  from  E    to  Eq1 j.   Typical  increases  in  excess  pore 
pressure  in  this  region  were  from  0.762  to  0.839,  from  0.568 
to  0.771;  from  0.591  to  0.688  and  from  0.445  to  0.614. 
Dissipation  then  proceeds  with  time  as  a  result  of  water 
flow  out  of  the  top  drainage  boundary,  Figure  5.36,  until 
time  factor  0.24632.   At  this  time  elements  48  and  53, 
adjacent  to  the  grid  center  line  in  the  third  row,  yield. 
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The  pre-  and  post-yield  pore  pressure  distributions  are  given 
in  Figure  5.37.   The  top  diagram  shows  that  in  the  time  between 
the  initial  yielding  and  just  prior  to  this  new  yielding  the 
pore  pressure  distribution  had  stabilized.   That  is,  all  local- 
ized areas  of  high  pore  pressure  had  disappeared  and  the  dis- 
tribution shows  the  smooth  contour  lines  of  the  no  yield  behav- 
ior.  In  the  lower  diagram  the  yielded  element,  number  53, 
exhibits  the  characteristic  pore  pressure  increase,  up  from 
0.3362  to  0.5736.   Dissipation  continues  until  T  =  0.26919, 
Figure  5.38,  when  a  small  effect  of  the  previous  yielding  still 
remained  and  when  elements  43  and  58  of  the  third  row  yielded. 
Two  more  elements  of  row  three  yielded  at  time  factor  0.36629 
again  demonstrating  the  localized  pore  pressure  increase  due 
to  the  yielding  of  a  particular  element  (Figure  5.39). 

Initially,  the  surface  settlements  are  equivalent  to  those 
obtained  from  the  no  yield  program.   These  displacements  shown 
in  Figure  5.35  are  due  to  distortion  within  the  clay  soil  and 
are  not  accompanied  by  any  volume  change.   Subsequent  settlements 
are  of  larger  magnitude  than  for  the  no  yield  case,  a  consequence 
of  the  increased  compressibility  of  the  yielded  elements.   This 
is  shown  in  Figure  5.40  where  surface  settlements  at  three  points 
are  plotted  versus  time  factor.   Values  are  normalized  by  divid- 
ing actual  settlements  by  the  ultimate  center-line  settlement 
for  the  yield  case.   For  the  case  of  no  yielding  heave  occurs 
initially  outside  the  loaded  area  and  then  decreases  with  time. 
For  the  yield  case  heave  increases  from  the  initial  value  due 
to  the  yielding  of  elements  under  the  load  then  decreases  with 
time . 
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5 . 6   Field  Studies 
In  recent  years  several  instrumented  embankments  have 
been  constructed  on  soft  clay  soils  and  described  in  the  soil 
mechanics  literature.   Bjerrum  (1972)  summarized  much  of 
this  work  with  regard  to  the  preconsolidation  effect.   His 
observations  were: 

1.  Almost  all  clay  soils  exhibit  such  an  effect 

2.  For  small  loads  within  the  range  of  the  effect 

a.  settlements  are  small 

b.  excess  pore  pressures  are  small 

c.  excess  pore  pressures  dissipate  rapidly. 

3.  For  larger  loads  which  exceed  the  range  of  the 
effect 

a.  settlements  are  large 

b.  excess  pore  pressures  are  large 

c.  excess  pore  pressures  dissipate  slowly. 

As  suggested  above,  and  in  most  reported  case  studies, 
the  nature  of  the  response  of  a  clay  foundation  to  a  super- 
imposed fill  is  related  to  the  weight  of  the  fill.   Figure 
5.41  is  reproduced  from  a  paper  by  Hoeg  et  al.  (1969)  and 
illustrates  the  variations  in  undrained  excess  pore  pressure 
at  different  locations  with  surface  pressure.   The  slope 
change  in  each  plot  indicates  yielding  of  the  soil  adjacent 
to  that  piezometer.   The  current  research,  however,  treats 
the  quasi-preconsolidation  phenomenon  more  fundamentally  in 
terms  of  effective  stresses,  and  illustrates  that  yielding 
may  occur  in  response  to  changes  in  both  isotropic  effective 
stresses  and  shear  stresses  independently  of  drainage 
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conditions.   When  a  load  is  applied  to  the  surface  of  the 
two-dimensional  problem  part  of  the  applied  total  stress  is 
carried  by  the  pore  water  and  part  by  the  soil  skeleton  as 
an  effective  stress.   If  the  applied  load  is  small,  the 
effective  stress  and  shear  stress  increments  will  be  small 
and  nowhere  will  the  effective  stress  state  exceed  the  yield 
surface.   If  however,  the  initial  load  is  larger,  or  if 
additional  shear  stresses  are  applied  the  stress  state  may 
lie  outside  the  yield  surface.   This  yielding  would  be 
expected  first  either  directly  under  the  load  and  subse- 
quently at  greater  depths  and  distances  from  the  load,  or 
at  a  point  near  the  edge  of  the  fill  where  shear  stresses 
are  high.   In  the  case  reported  by  Hoeg  et  al.  point  A  on 
the  fill  center  line  yielded  first,  followed  at  greater 
loads  by  points  B  and  C.   Yielding  under  an  initial  load 
was  illustrated  in  the  finite  element  two-dimensional  study 
in  Figure  5.35. 

An  insitu  yield  curve  may  be  back  figured  for  the  partic- 
ular depth  of  the  piezometers  in  Figure  5.41.   From  the  known 
surface  loading,  total  stresses  and  consequently  effective 
stresses  can  be  calculated  at  the  piezometer  locations. 

Total  stress  increments  were  calculated  assuming  a  uni- 
form circular  load  on  a  homogeneous  elastic  half  space. 
Tabulated  stress  functions  were  obtained  from  Ahlvin  and 
Ulery  (1962).   In  Figure  5.42  stress  paths  for  the  four 
locations  are  plotted  in  terms  of  the  effective  octahedral 
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stress  increments.   The  yield  points  define  a  yield  curve 
similar  in  shape  to  that  found  in  the  laboratory. 

Yielding  may  also  occur  during  consolidation.   As  excess 
pore  pressures  decrease  due  to  flow  of  water  out  of  the  clay 
stratum,  effective  stresses  increase  and  may  at  a  point 
exceed  the  yield  surface.   This  effect  was  observed  in  the 
one-,  two-,  seventy  five-,  and  one  hundred-element  studies. 
These  solutions  also  illustrate  that  the  foundation  yielding 
is  progressive  rather  than  occurring  at  one  particular  load. 
Zones  of  the  soil  may  yield  due  to  initial  or  subsequent 
loadings,  while  other  zones  may  yield  progressively  during 
consolidation. 

Results  from  the  finite  element  study  may  be  summarized 
as  follows.   If  the  incremental  effective  stress  state  for 
an  applied  load  plots  within  the  yield  surface,  then  the 
settlements  are  small,  pore  pressures  are  small  and  dissipate 
rapidly.   Should  the  effective  stress  state  cross  the  yield 
surface,  due  to  additional  loading  or  consolidation,  then 
settlements  become  large,  pore  pressure  increases  due  to 
load  transfer  from  skeleton  to  pore  water,  and  dissipation 
is  slow. 

The  overall  soil  behavior  as  predicted  by  the  finite 
element  study  is  in  general  agreement  with  reported  field 
observations . 
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CONCLUSIONS  AND  RECOMMENDATIONS 


6.1   Conclusions 


In  this  dissertation  laboratory  and  analytical  pro- 
cedures were  combined  to  study  the  effects  of  quasi  pre- 
consolidation  on  the  compressibility  of  clay  soils.   A  soil 
'behavior  model  was  formulated  from  the  results  of  triaxial 
tests  performed  on  clay  specimens  in  which  a  quasi  precon- 
solidation  effect  had  been  induced.   The  model  was  then 
introduced  into  the  finite  element  procedure  such  that 
boundary  value  consolidation  problems  could  be  solved.   From 
the  study  performed  the  following  conclusions  may  be  drawn. 

1.  A  quasi  preconsolidation  effect  can  be  induced  in 
a  clay  specimen  by  consolidating  it  in  a  triaxial  cell, 
under  a  hydrostatic  confining  pressure.   Loading  periods 
of  approximately  21  days  were  employed.   This  procedure 
eliminates  the  disturbance  due  to  stress  release  and  sample 
preparation  experienced  when  samples  are  trimmed  from  blocks 
prepared  in  a  consolidometer  type  apparatus. 

2.  Drained  triaxial  tests  performed  on  such  specimens 
along  different  stress  paths,  show  essentially  similar 
stress-strain  behavior.   In  every  case  the  plot  is  initially 
steep,  followed  by  a  fairly  sharp  transition  curve  to  an 
approximately  straight  line  with  shallow  slope  At  this 
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transition,  the  compressibility  increases  abruptly. 

3.  From  the  triaxial  test  results  a  soil  model  may  be 
formulated  with  the  following  features. 

a.  The  stress-strain  behavior  can  be  represented  by 
two  straight  lines  intersecting  at  a  point,  termed 
the  q-yicld  point. 

b.  For  all  tests  performed  on  specimens  prepared  under 
the  same  confinement  pressure,  the  q-yield  points 
define  a  curve  which  can  be  closely  approximated 

by  an  ellipse. 

c.  The  yield  point  ellipses  for  different  confinement 
pressures  are  found  to  form  cross  sections  of  an 
elliptic  cone  with  apex  at  the  zero  confinement 
origin. 

This  three-dimensional  yield  surface  can  be  defined  by 
a  single  equation. 

4.  The  soil  model  derived  from  laboratory  test  results, 
for  a  clay  soil  which  exhibits  a  quasi  preconsolidation 
effect,  can  be  introduced  into  the  finite  element  procedure. 
Due  to  the  fundamental  versatility  of  the  finite  element 
procedure,  consolidation  solutions  can  be  obtained  for  a 
great  variety  of  plane  strain  configurations. 

5.  Solutions  were  obtained  for  one-  and  two-dimensional 
problems  with  and  without  consideration  of  the  quasi  precon- 
solidation behavior  for  a  number  of  finite  element  mesh 
sizes.   Results  illustrate 
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a.   fundamental]  behavior  at  a  point  within  a  clay  soil 
prior  to,  ;it,  and  after  the  quasi  preconsolidation 
pressure  i :;  exceeded.   When  yielding  occurs  "the 
compressibility  of  the  soil  skeleton  increases  and 
a  portion  of  the  stress  carried  by  it  as  an  effective 
stress  is  transferred  to  the  pore  water.   Subsequent 
settlements  are  increased  and  the  rate  of  pore 
pressure  dissipation  decreased  relative  to  no 
yielding  behavior, 
and  b.   the  overall  effects  of  construction  on  a  clay 

stratum  which  exhibits  the  quasi  preconsolidation 
effect.   Effective  stresses  from  construction 
loads  which  do  not  exceed  the  insitu  effective 
stress  yield  surface  result  in  small  settlements 
and  rapid  dissipation  of  excess  pore  pressures. 
If  yielding  does  occur  then  the  corresponding 
settlements  are  large  and  excess  pore  pressure 
dissipation  slow.   These  are  in  agreement  with 
field  observations  reported  in  the  literature. 

6. 2   Recommendations 
In  the  current  research  the  laboratory  study  was  per- 
formed to  establish  a  utilitarian  soil  model  that  might  be 
incorporated  into  the  finite  element  procedure.   For  a  more 
definitive  investigation  of  the  quasi  preconsolidation 
effect  under  a  three  dimensional  stress  state  the  following 
recommendations  are  made. 
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1.  Further  experiment  using  commercial  soils,  such  as 
grundite,  should  be  conducted: 

a.  to  investigate  the  effects  of  length  of  q-con- 
finement  time 

b.  to  establish  the  uniqueness  of  the  yield  curve 
and  its  shape  for  negative  Aa   .  . 

c.  to  verify  the  linear  cone  relationship  by 
considering  a  greater  range  in  confinement 
pressures,  especially  at  low  values. 

2.  The  model  for  the  quasi-preconsolidation  effect 
should  be  verified  for  specimens  consolidated  under 
KQ  conditions.   Comparison  should  be  made  with 
tests  performed  on  specimens  trimmed  from  clay 
consolidated  in  an  oedometer  type  apparatus. 
Effects  of  anisotropy  should  be  considered. 

3.  The  applicability  of  the  model  to  clays  consolidated 
in  the  field  should  be  studied  using  high  quality 
undisturbed  samples. 

The  model  established  for  the  clay  soil  in  this  research 
was  that  of  a  bi-linear  elastic  material,  with  the  inter- 
section point  defined  by  a  three-dimensional  effective  stress 
surface.   The  finite  element  procedure  is  of  such  versatility 
that  more  complicated  behavior  corresponding  more  closely  to 
real  response,  can  be  accommodated.   The  soil  behavior  might  be 
more  accurately  modeled  using  incremental  stress-strain  curves, 
by  considering  anisotropic  conditions  and  by  including  an 
ultimate  failure  criterion. 
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The  actual  finite  element  procedure  used,  proved  to  be 
extremely  costly  in  computer  time,  due  to  the  exceedingly 
short  time  increments  which  must  be  employed.   Programs 
consisting  of  100  elements  achieved  approximately  200  cycles 
of  dissipation  in  1,000  seconds  of  CDC  computer  time.   When 
determining  stresses  and  deformations  for  field  situations 
using  the  finite  element  method,  600  to  1000  elements  are 
typically  employed.   Such  a  grid  for  the  consolidation  prob- 
lem would  result  in  exorbitant  computer  running  times.   A 
need  exists,  therefore  for  a  more  economical  finite  element 
procedure.   More  efficient  storage  to  reduce  memory  require- 
ments and  the  possibility  of  an  implicit  rather  than  explicit 
solution  procedure  might  be  considered. 

Prediction-observation  studies  of  field  structures  are 
strongly  recommended.   Solutions  from  computer  programs  includ- 
ing models  based  on  laboratory  testing  of  field  samples  need 
to  be  compared  with  field  instrumentation  measurements. 

Finally,  the  question  arises  as  to  the  influence  of  the 
quasi-preconsolidation  pressure  on  strength  and  how  strength 
parameters  should  be  determined  in  the  laboratory.   The 
implications  of  this  question  to  stability  analyses  should 
be  investigated. 
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A.   SAMPLE  PREPARATION 

A.  1   Sedimentation  Apparatus 

A  steel  frame  designed  and  constructed  at  Purdue  Univer- 
sity by  Altschaeffl  (1960)  was  used  to  sediment  the  samples. 
A  general  view  of  the  eight  bay  frame  is  shown  in  Figure  A.l 
•and  a  line  diagram  indicating  the  arrangement  in  one  bay  is 
shown  in  Figure  A. 2.   The  front  tank  was  filled  with  petrola- 
tum oil  (specific  gravity  0.84)  such  that  it  balanced  exactly 
the  piston  and  hanger  with  four  16  kg.  counter  weights  in 
place.   The  connecting  wire  cable  was  1/8  inch  diameter, 
2000  lb.  test  cable  passing  over  balanced  twelve-inch  diame- 
ter ball  bearing  pulleys. 

The  lucite  tube  in  which  the  sample  was  to  be  sedimented 
was  20  inches  long  and  3-1/2  inches  in  diameter.   A  6  inch 
diameter  by  1  inch  thick  lucite  base  with  drainage  outlet  was 
screwed  to  the  flange  on  the  sedimentation  tube.   The  loading 
piston  was  also  made  of  lucite  and  had  a  porous  stone  and 
alignment  spiders  attached.   A  close  up  of  the  piston-tube 
apparatus  is  shown  in  Figure  A. 3. 

Loading  of  the  sample  was  carried  out  by  transferring 
oil  by  means  of  a  capillary  siphon  from  the  front  tank  into 
a  second  completely  independent  rear  tank.   The  load  on  the 
sample  was  determined  by  lowering  the  rear  tank  onto  a  100  kg 
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FIGURE    A- 1     THE     SEDIMENTATION     FRAME    . 
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FIGURE     A-3    LUCITE     TUBE     AND     PISTON. 
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capacity  platform  scale  (sensitivity  10  gms . )  and  weighing 
the  oil  transferred.   Capillary  siphons  were  made  from  glass 
tubing  with  bore  diameters  ranging  from  0.75  to  3  mms . 
Dimensions  are  shown  in  Figure  A. 4.   All  siphons  were  ini- 
tially numbered  and  calibrated  with  a  12  inch  unbalanced 
head  (Table  A.l) . 

A. 2   Sedimentation  Procedure 
The  Grundite  was  first  dry  sieved  and  only  that  portion 
which  passed  through  a  U.  S.  Bureau  of  Standards  Number  200 
sieve,  opening  0.074  mm,  used.   A  Calgon  solution,  to  be 
used  as  a  clay  particle  dispersant,  was  prepared  by  dissolv- 
ing 41  gms  of  commercial  Calgon  in  1000  ccs  of  distilled 
water.   Three  batches  of  slurry  were  made,  each  consisting 
of  250  gms  of  Grundite,  125  ccs  of  Calgon  solution  and  375 
ccs  of  distilled  water.   These  were  mixed  in  milkshake 
machines  for  half  an  hour.   Six  drops  of  0.1  N  IIC1 ,  a  parti- 
cle flocculant,  were  added  to  each  batch  one  minute  before 
mixing  terminated.   The  slurries  were  then  washed  into  a 
large  glass  jar  using  an  additional  500  ccs  of  distilled 
water.   The  preparation  apparatus  is  shown  in  Figure  A. 5. 
The  jar  was  attached  to  a  vacuum  and  the  mixture  dcaircd  and 
stirred  for  a  further  half  hour.   The  jar  was  then  inverted 
and  the  slurry,  consisting  of  750  gms  of  Grundite  in  2000 
ccs  of  water,  tremied  into  the  lucite  sedimentation  tube. 
This  procedure  was  used  so  as  to  minimize  the  possibility  of 
trapping  air  bubbles  in  the  sample. 
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FIGURE      A-4      GLASS      CAPILLARY    SIPHON 
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TABLE  A.l 
Glass  Siphon  Calibrations 

Grams  of  oil  transferred/24  hours 

Siphon  Number  (12  inch  unbalanced  head) 

1  990 

2  1170 

3  1090 

4  980 

5  680 

6  810 

7  400 

8  1050 

9  1000 

10  1150 

11  140 

12  250 

13  260 

14  310 

15  320 

16  4000 

17  2500 

18  2810 

19  2750 
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FIGURE    A-5     PREPARATION     EQUIPMENT    FOR    SEDIMENTATION 
SAMPLES. 
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The  tube  had  been  prepared  as  follows.   First  the  in- 
side walls  were  sprayed  with  a  fine  coat  of  Silicone  using 
a  Scars  Silicone  Spray  aerosol.   The  tube  flange  and  base 
rim  were  greased  to  prevent  leakage.   The  porous  stone, 
carefully  deaired  and  covered  with  filter  paper,  was  placed 
on  the  base  and  the  tube  screwed  into  place.   After  the 
slurry  had  been  tremied  into  the  lucite  tube  the  soil  was 
allowed  to  settle  for  48  hours.   The  drainage  outlet  in  the 
base  was  kept  open  to  create  a  small  downward  flow  of  water 
which  caused  some  consolidation  of  the  clay.   Figure  A. 6 
shows  a  freshly  prepared  sample  and  one  after  the  48  hour 
standing  period.   In  early  samples  consolidation  was  further 
assisted  by  applying  a  vacuum  at  the  drainage  outlet.   This 
was  discontinued  as  satisfactory  results  were  obtained 
during  the  settling  period  alone. 

The  lucite  tube  was  then  placed  in  the  sedimentation 
frame,  the  piston  porous  stone  deaired  and  placed  in  contact 
with  the  soil,  and  the  oil  tank  balanced.   A  capillary 
siphon  with  very  slow  rate  of  loading  was  then  introduced 
and  the  rear  tank  weighed  every  24  hours.   As  the  soil  con- 
solidated and  its  strength  increased  it  was  possible,  after 
a  few  days,  to  substitute  a  larger  bore  siphon  or  add  an 
additional  siphon.   After  sufficient  strength  gain  incre- 
mental loads  of  4  kg.  were  added  to  the  piston  hanger.   A 
typical  loading  plot  is  shown  in  Figure  A. 7.   Samples  were 
all  loaded  to  50  kg.  (11.45  psi)  in  about  three  weeks.   They 
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FIGURE     A-6     SEDIMENTATION     SAMPLES    BEFORE    AND 
AFTER    STANDING    PERIOD   . 
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remained  under  this  load  for  three  days  before  being  ex- 
truded for  triaxial  sampling. 

Homogeneity  of  each  sample  was  checked  in  an  approxi- 
mate manner  by  determining  the  water  content  of  specimens 
from  the  upper,  lower  and  middle  third  regions  of  the 
extruded  sample.   In  no  case  was  the  moisture  content  varia- 
tion more  than  10%.   Hydrometer  analyses  were  also  performed 
on  three  specimens  from  one  sample  and  the  grain  size  dis- 
tribution plots  are  given  in  Figure  A. 8.   No  significant 
difference  in  composition  is  apparent.   Every  effort  was 
made  to  standardize  preparation  procedures  so  as  to  minimize 
sample  variation.   Void  ratio  determinations  were  made  for 
all  samples  and  a  range  of  values  from  1.52  to  1.62  with  an 
average  of  1.57  was  found. 

A. 3   Triaxial  Apparatus 
A. 3.1   The  Triaxial  Cell  Set-up .   In  order  to  maintain 
a  constant  confining  pressure  on  triaxial  samples  for  a  long 
period  of  time  the  following  set-up  was  designed.   The  sys- 
tem is  shown  in  Figure  A. 9  and  a  line  diagram  showing  the 
essential  arrangement  is  given  in  Figure  A. 10.   The  air 
pressure  was  from  a  Curtis  two-stage  compressor  with  a  100 
psi  capacity.   A  back-up  Nitrogen  tank  was  available  via  a 
ball  valve  should  the  primary  air  pressure  drop.   Three 
Fairchild  Miller  multi  stage  pressure  regulators  and  dials 
called  SOURCE  I,  SOURCE  II  and  SOURCE  III  were  connected  in 
parallel  with  the  capability  of  inserting  further  portable 
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SEDIMENTATION    SAMPLE. 
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sources  if  required.   A  loop  of  1/8  inch  inner  diameter 
plastic  tubing  had  connections  into  which  eight  triaxial 
cells  could  be  attached.   The  pressure  in  the  loop  might  be 
supplied  by  one  source  or  the  loop  could  have  been  broken 
at  any  one  of  several  quick  connectors  and  different  pres- 
sures applied  to  different  cells.   When  a  triaxial  sample 
was  to  be  tested  a  length  of  tubing  with  connectors  at  each 
end  would  be  inserted  and  the  sample  tested  in  another  area 
of  the  room  while  still  under  the  same  confining  pressure. 

A. 3. 2   The  Triaxial  Cells .   The  following  makes  of 
triaxial  cells  were  used  in  the  experimentation:  - 

1.  Geonor  cell,  Norwegian  Geotechnical  Institute 

2.  Wykeham  cell,  Wykeham  Farrance  Eng .  Ltd.,  Slough, 
England 

3.  Test  Lab  cell,  Test  Lab.  Corporation,  Chicago, 
Illinois 

The  Geonor  apparatus  has  been  described  in  detail  by 
Andresen  and  Simons  (1960)  and  was  the  cell  used  primarily 
in  this  research.   Essentially  it  consists  of  two  main 
parts,  a  base  plate  with  pedestal  on  which  the  sample  sits 
and  an  upper  plate  with  surrounding  plexiglass  cylinder. 
These  parts  are  screwed  together  by  three  long  wing  bolts, 
as  shown  in  Figure  A. 11.   Two  outlets  in  the  base  pedestal 
permit  drainage  from  the  bottom  of  the  sample  to  occur  or 
can  be  used  for  the  measurement  of  the  pore  water  pressure. 
Drainage  from  the  top  of  the  sample  is  possible  by  connect- 
ing a  special  platen  to  another  outlet  in  the  cell  base. 
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FIGURE     A-ll      THE     GENOR    TRIAXIAL    TEST. 


FIGURE     A-12     GEONOR     CELL ,  ROTATING    BUSHING     DETAIL. 
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The  Geonor  cells  used  were  all  equipped  with  rotating 
bushings  to  reduce  piston  friction.   Figure  A. 12  shows  the 
assembly,  and  Figure  A. 13  a  plot  from  Andrescn  and  Simons 
(1960)  indicating  the  considerable  advantage  of  using  the 
rotating  bushing.   Simple  experiments  were  run  "testing"  a 
stiff  rubber  block  and  a  spring  as  shown  in  Figure  A. 14. 
Results  were  not  as  spectacular  as  those  given  by  Andresen 
and  Simons  but  tests  run  without  the  rotating  bushing  did 
show  occasional  jagged  breaks  on  the  output  curve  due  to 
piston  sticking  and  releasing. 

The  Geonor  cell  was  also  adapted  for  use  in  a  series  of 
isotropic  loading  tests.   In  these  tests  the  cell  confining 
pressure  was  increased  incrementally  at  regular  time  inter- 
vals.  It  was  desired  to  measure  both  axial  deformation  and 
volumetric  change.   The  latter  was  easily  accomplished  by 
permitting  drainage  from  both  ends  of  the  sample  and  attach- 
ing the  outlet  lines  to  a  calibrated  burette.   The  axial 
deformation  was  measured  in  the  following  manner.   A  linear 
variable  differential  transformer  (LVDT)  core  was  threaded 
and  screwed  into  a  top  platen.   The  LVDT  itself  was  attached 
to  the  piston  which  was  fixed  to  the  cell  and  not  permitted 
to  move.   Any  axial  movement  of  the  top  platen  produced  a 
corresponding  movement  of  the  core  within  the  fixed  LVDT  and 
the  output  was  recorded  automatically.   Figure  A. 15  shows  a 
photograph  of  the  set-up,  in  this  case  with  only  bottom  end 
drainage  permitted.   The  leads  from  the  LVDT  were  coated 
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FIGURE     A-15     SET-UP    FOR    ISOTROPIC     LOADING    TEST. 


17  7 

with  plastic  and  exited  through  a  carefully  sealed  hole  cut 
in  the  cell  base.   Calibration  was  accomplished  by  moving 
the  piston  and  sighting  a  point  on  the  LVDT  with  a  cathetome 
ter.   These  actual  measurements  were  related  to  the  number 
of  divisions  transcribed  on  the  recorder  output. 

The  Wykeham-Farrance  cell  is  shown  in  Figure  A. 16. 
This  apparatus,  like  the  Geonor,  consists  of  a  base  plate 
with  pedestal  and  an  upper  plate  with  plexiglass  cylinder. 
Double  drainage  of  the  soil  sample  is  also  permitted.   The 
main  disadvantage  of  this  cell  is  the  absence  of  a  rotating 
bushing;  consequently  piston  friction  is  an  unknown  factor. 

A  few  tests  were  carried  out  in  a  Test  Lab  triaxial 
cell.   This  apparatus  shown  in  Figure  A. 17  is  very  similar 
to  the  Wykeham  Farrance  cell  and  has  the  same  frictional 
disadvantages . 

A. 3. 3   Testing  Machine  and  Recorders .   The  testing 
machine  used  for  all  the  constant  strain  rate  tests,  and 
shown  in  Figure  A. 18,  was  manufactured  by  Engineering  Labora- 
tory Equipment  Ltd.  (ELE)  of  Rickmansworth ,  England.   The 
soil  sample  within  the  triaxial  apparatus  was  compressed  by 
raising  the  platform  on  which  the  cell  sits  while  the  piston 
bore  against  a  proving  ring  or  load-cell.   In  this  research 
a  100  lb.  capacity  Baldwin-Lima-Hamilton  Corporation  load- 
cell  was  used.   The  E.L.E.  machine  allows  a  wide  range  of 
testing  rates  varying  from  1.6  x  10"1  to  2.0  x  10"e  inches 
per  minute.   For  a  2.8  inch  sample  this  represents  a  range 
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FIGURE     A-16     THE     WYKEHAM-FARRANCE     CELL 


FIGURE    A-17    THE     TEST    LAB     CELL 
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FIGURE     A-18    THE     TRIAXIAL    TESTING    MACHINE 
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of  one  percent  strain  in  10.5  seconds  to  one  percent  strain 
in  approximately  10  days. 

The  recorder  used  during  triaxial  testing  for  both  load- 
cell  and  volume  change  output  was  a  Sanborn  Carrier  Amplifier 
Recorder,  Model  321.   This  is  a  portable,  two-channel  record- 
ing system  using  hot-wire  writing  arms  on  heat  sensitive 
paper.   It  is  shown  in  Figure  A. 19.   For  the  Linear  Variable 
Differential  Transformer  output  in  the  Isotropic  loading 
tests  a  single  channel  ink  recorder,  Varian  Model  G-14A-1 
was  employed  (Figure  A. 20). 


APPENDIX  B 
LABORATORY  TEST  RESULTS 
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FIGURE     A-19     THE    SANDBORN     321    RECORDER 


FIGURE     A-20    THE     VARIAN     G-I4A-I     RECORDER 
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B.   LABORATORY  TEST  RESULTS 


B.l   Introduction 


A  particular  test,  Number  23,  is  used  as  an  example  and 
is  described  in  detail  in  the  following  section.   Results 
from  all  phases  of  its  preparation  and  testing  are  given. 
•Section  B.3  contains  detailed  results  from  all  twenty-seven 
tests . 

B.2   Triaxial  Test  Number  23 
Sedimentation  sample  ID  was  prepared  on  June  7,  1972  as 
described  in  Section  A. 2  and  allowed  to  consolidate  under 
its  own  weight  with  the  base  drainage  outlet  open.   Forty- 
eight  hours  later  the  lucite  tube  was  placed  in  the  sedimen- 
tation frame  and  the  piston  and  deaired  porous  stone  brought 
into  contact  with  the  soil  surface.   A  reference  point  was 
marked  on  the  piston  such  that  cathetometer  readings  of  the 
sample  compression  could  be  measured  during  loading.   Siphon 
number  11  was  inserted  and  rear  tank  weight  and  cathetometer 
readings  taken  every  24  hours.   After  twelve  days  siphon 
number  7  with  a  slightly  larger  bore  was  substituted.   Still 
later  siphons  9  and  17  were  used  and  finally  4  kg.  load 
increments  were  added.   The  final  load  of  50  kg.  is  equiva- 
lent to  a  vertical  pressure  of  11.4  psi.   Figure  B.l  shows  a 
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plot  of  the  loading  versus  time,  and  Figure  B.2  the  e  -  log  p 
curve  obtained.   The  final  load  increment  was  added  on  July 
10th  and  allowed  to  remain  for  three  days. 

On  July  13th  the  soil  was  extruded,  six  water  contents 
of  58,  58,  59,  59,  59,  and  60  obtained  at  different  depths 
within  the  clay  and  a  specimen  cut  for  Triaxial  Test  Number 
23.   The  summary  sheet  for  this  test  is  shown  in  Figure  B.3. 
The  triaxial  cell  was  placed  in  the  triaxial  set-up  in  posi- 
tion VII  C  at  a  confining  pressure  of  30  psi.   The  drainage 
lines  were  connected  to  a  burette  and  volume  change  readings 
taken  at  intervals  during  the  q-confining  period.   These 
results  are  shown  in  Table  B.l.   From  these  values,  the 
excess  water  (2.10  cc.)  contained  initially  in  the  filter 
paper  which  encircled  the  specimen,  must  be  subtracted.   On 
August  4,  1972  when  the  sample  was  to  be  tested  the  total 
volume  change  which  had  occurred  during  the  twenty- two  day 
confining  period  was  13.00  cc.   Since  there  was  no  discernible 
visual  distortion  of  the  sample,  new  dimensions  were  calculated 
assuming  that  the  sample  shape  remained  a  right  cylinder  with 
the  same  length-diameter  ratio. 

ttD2 
Sample  volume     Vol  =  —j—  L 


Assumption   r  =  D/L  =  constant 

Vol  =  £J-  L3  =  0.1830  L3 
Original  sample  volume  =  73.15  ccs 
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FIGURE     B-2      VOID    RATIO -PRESSURE    CURVE    FOR 
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Test    Number 23.... 

Quasi  Number VJI.C. 


John  L  Davidson 
101  H  Grissom   Hall 


Sedimentation    Sample  Prepared  .31:72. 

Loading   Started  ...6,9,72 

Quasi  Sample   Set   Up       .7,13,72 

Quasi   Pressure   Applied   ..?.-. 1 3.- 72. 

Triaxial  Test   Started  8-.4-.72 


Initial    Void    Ratio   ...157.... 

Unit    Weight. ..JQ8 pcf. 


r  I  Top   42.1 
2        42.2 


Water    Content 
after   Test 


Volume   Change    During   Quasi   Set-Up !3,00 . ....ccs 

Prior  to  Testing  Sample    Length.... 2.72 ins. 

Sample  Area... .1.351.0. sq.ins. 


3 

4 

5 

"-  6 


42.1 
42.3 
42.2 
42.2 


Cell  Type  ...WYKENAM 


Standard    Procedures  7   ..YES. 


At. 


&V, 


oct 


5.92 


psi 


^T-oct     c  -  -8^-6  ---  psi 


Stress    Path 


FIGURE    B-3    SAMPLE    SUMMARY    SHEET. 
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TABLE  B.l 

Volume  Change  during  the  q-consolidation  Period 
Test  Number  23 


Date 

Burette 

AV 

7 

13 

72 

1.30 

0.00 

7 

14 

72 

9.10  -  0.18 

7.80 

7 

15 

72 

3.00 

10.62 

7 

16 

72 

3.80 

11.42 

7 

17 

72 

5.30 

12.92 

7 

18 

72 

5.36 

12.98 

7 

19 

72 

6.40 

14.02 

7 

20 

72 

2.18  -  6.53 

14.15 

7 

21 

72 

2.26 

14.23 

7 

24 

72 

2.52 

14.49 

7 

26 

72 

2.66 

14.63 

7 

28 

72 

2.75 

14.72 

7 

31 

72 

2.90 

14.87 

8 

2 

72 

3.05 

15.02 

8 

4 

72 

3.13 

15.10 
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New  sample  volume  =  73.15  -  13.00  =  60.15  ccs 

=  0.1830  L"  ■ 
new 

L  ,,  =  6.90  cms 
new 

=2.72  inches 

D    =1.31  inches 
new 

Area    =  1.3510  sq  ins 
new  n 

These  calculations  were  carried  out  only  in  Triaxial  Test 
Number  23.   Consequently,  in  all  other  tests  there  is  an 
"underprediction  in  axial  strain,  axial  stress  and  volumetric 
strain  due  to  use  of  larger  than  actual  sample  dimensions. 

Triaxial  Test  Number  23  was  a  Type  1  test  with  a  constant 
chamber  pressure  of  30  psi,  equal  to  the  q-consolidation 
pressure.   An  extension  plastic  tube  with  quick  connectors 
at  each  end  was  inserted  so  that  the  sample  could  be  placed 
in  the  loading  frame  and  still  be  attached  to  position  VII  C 
at  the  q-consolidation  pressure.   The  burette  was  removed  and 
the  drainage  lines  connected  to  the  special  burette  with 
transducer  attachment  described  in  Section  3.6.   The  initial 
test  parameters  were: 

Triaxial  machine  loading  rate:  4x10" 5  inches/minute 

Sandborn  chart  speed:    1.2  divisions/minute 

Load  cell  calibration:    0.4  lb. /division 

Volume  change  burette  calibration:  ~  0.05  cc. /division 

It  was  found  to  be  more  convenient  to  calibrate  the  burette 
accurately  during  the  test.   Burette  levels  were  read 
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visually  at  different  times  throughout  the  test  and  recorded 
directly  on  the  Sandborn  chart.   After  the  test  these  read- 
ings were  compared  with  the  transducer  output  and  the  best 
calibration  constant  derived.   For  Test  23  this  was  0.0464 
cc. /division .   During  the  test  the  Sandborn  attenuator  was 
altered  to  reduce  the  output  scale,  changing  the  load  cell 
calibration  to  0.8  lb ./division .   The  test  was  continued  at 
constant  strain  rate  until  approximately  15  percent  axial 
strain  had  been  achieved.   Finally  the  triaxial  sample  was 
removed  and  water  content  determinations  made. 

Test  data  consisted  of  the  Sandborn  chart  with  plots  of 
load  cell  output  and  volume  change  burette  output.   A  computer 
program  was  written  to  facilitate  the  reduction  of  the  test 
data.   After  reading  in  the  initial  test  parameters  described 
above,  it  was  only  necessary  to  punch  cards  containing  three 
variables.   These  were  the  number  of  divisions  along  the 
chart,  number  of  load  cell  divisions  and  number  of  burette 
divisions.   Computer  output  consisted  of  listings  of  axial 
strain,  axial  stress,  volumetric  strain,  deviatoric  strain, 
octahedral  normal  and  shear  stresses,  axial  load,  axial 
deformation  and  corrected  sample  areas.   The  Calcomp  Plotter 
was  then  used  to  plot  the  desired  stress-strain  relationships. 
A  program  description,  listing,  flow  chart  and  example  prob- 
lem are  given  in  Appendix  C.   Figures  (B.4),  (B.5)  and  (B.6) 

show  respectively  the  At   .  vs  e,  At   .  vs  e   .  , ,  and 
r       '  oct      '    oct     axial' 

v  vs  Aa    plots  for  Test  Number  23.   These  figures  were  traced 
from  the  Calcomp  plots  to  provide  visual  uniformity. 
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B.3  Triaxial  Test  Plots 
Figures  B.4  through  B.41  give  results  from  the  twenty- 
seven  triaxial  tests  performed  in  the  laboratory  study. 
Stress  and  strain  parameters  were  defined  in  Section  3.5. 
A  brief  summary  of  the  parameters  used  as  axes  descriptors 
in  the  following  figures  is  now  given.   Note  that  all  stresses 
are  increments  from  the  qc_  (octahedral  normal  stress  prior 
to  shearing)  values  given  in  Table  4.2. 

Axial  deformation  -  change  in  specimen  length  [  AL  ] 

Axial  strain  -  change  in  specimen  length  divided  by  original 
length  [  AL/L  ] 

Axial  stress  -  axial  load  divided  by  specimen  cross  sec- 
tional area  [  P/As  ] 

Cell  pressure  -  triaxial  cell  fluid  pressure  [  a   ] 

Deviatoric  strain  -  major  principal  strain  minus  one-third 
the  volumetric  strain  [  Acj  -  6v/3  ] 

Deviator  stress  -  difference  between  major  and  minor  prin- 
cipal stresses  [  AOj  -  Aa3  ] 

Octahedral  normal  stress  -  mean  of  principal  effective 
stresses  [  (Ao1  +  2Aa3)/3  ] 

Octahedral  shear  stress  -  factor  times  principal  effective 
stress  differences  [  /2(Aaj  -  Aa3)/3  ] 

Volume  change  -  specimen  volume  change  [  AV  ] 
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Volumetric  strain  -  volume  change  divided  by  original 
volume  [  AV/V  ] 
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.o  I  1.5 

Volume     Change    (cc) 


Axial    Deformation    (xlO     in) 
FIGURE    B-22    TEST    NUMBER    14    AXIAL    STRESS    INCREMENT- 
VOLUME    CHANGE    AND    AXIAL    DEFORMATION  . 
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10  15         „  20 

Axial    Deformation    (x!0"^in) 
FIGURE     B-23    TEST    NUMBER    15    OCTAHEDRAL    SHEAR 
STRESS  -  VOLUME     CHANGE     AND     AXIAL 
DEFORMATIONS      PLOTS. 
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0.4  0.6  0.8 

Volume    Change    (cc) 


.2 


4  6 

Axial    Deformation 


(x  I0"4  in) 


FIGURE     B-26 


TEST    NUMBER    IS    OCTAHEDRAL    SHEAR 
STRESS  -  VOLUME    CHANGE    AND    AXIAL 
DEFORMATION    PLOTS. 
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6  8 

Axial     Deformotion 


FIGURE     B-34    TEST    NUMBER    22    CELL    PRESSURE- 
AXIAL    DEFORMATION     PLOT  . 
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THE  TRIAXIAL  TEST  REDUCTION  PROGRAM 


C.l   Introduction 


Load  cell  and  sample  volume  data  from  drained  triaxial 
tests  are  recorded  as  continuous  curves  on  a  single  roll  of 
output  paper  by  the  Sandborn  recorder.   This  computer  pro- 
gram has  been  written  to  facilitate  the  meaningful  interpre- 
tation of  the  output  by  calculating  relevant  stresses  and 
strains  and  listing  these  in  tabular  form  and  if  desired  in 
CALCOMP  plots. 

A  flow  chart  for  the  program  is  given  in  Figure  C.l. 
Section  C.2  contains  the  program  input  format  and  computer 
word  definitions.   Section  C.3  shows  a  simple  example  con- 
sisting of  only  thirteen  points.   Input,  output  and  plots 
are  shown.   A  listing  of  the  program  is  given  in  Section 
C.4.   Comment  cards  have  been  used  to  explain  the  step  by 
step  manipulations. 
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(start) 
"T" 


\  READ        ~~L 

Yest  FACTORS/ 

V"       READ        ~~7 
\  CHART    INPUT /- 


CALCULATE- 

LOAD 

DISPLACEMENT 
STRAIN 

CORRECTED    AREA 
CORRECTED   VOLUME 
AXIAL    STRESS 
VOLUME    CHANGE 
VOLUMETRIC   STRAIN 
DEVIATORIC    STRAIN 
OCTAHEDRAL  STRESSES 


NO 


YES 


FIGURE     C-l      FLOW    CHART    FOR    TRI AXIAL   TEST 
REDUCTION    PROGRAM. 
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C . 2   Input  Format  for  Triaxial  Program 

IDENTIFICATION  CARD  (3A10) 
Columns  1-30   TITLE  (3) 


ii.   PLOT  CARD  (18) 

Column  1    I  PLOT  =  0 
=  1 


Column 

2 

IPLOT 

= 

0 
2 

Column 

3 

I  PLOT 

= 

0 

3 

Column 

4 

IPLOT 

= 

0 

Column 

5 

IPLOT 

= 

0 
5 

Column 

6 

IPLOT 

= 

0 
6 

Column 

7 

IPLOT 

■ 

0 
7 

Column 

8 

IPLOT 

= 

0 

if  no  plot  required 

Oct.  Shear  Stress  -  Oct.  Normal 

Stress 
if  no  plot  required 
Oct.  Shear  Stress  -  Deviatoric 

Strain 
if  no  plot  required 
Oct.  Shear  Stress  -  Volumetric 

Strain 
if  no  plot  required 
Oct.  Shear  Stress  -  Axial  Strain 
if  no  plot  required 
Oct.  Normal  Stress  -  Axial  Strain 
if  no  plot  required 
Volumetric  Strain  -  Oct.  Normal 

Stress 
if  no  plot  required 
Deviator  Stress  -  Axial  Strain 
if  no  plot  required 
Volumetric  Strain  -  Deviatoric 

Strain 


in 


IV, 


TEST  FACTORS  (4F10.0) 

Columns   1-10    PC  -  Sandborn  chart  speed   divs/minute 
11-20    FL  -  Load-cell  calibration 

factor  lbs/diA^isi  on 

21-30         FT    =   Triaxial   machine 

loading  rate  ins/minute 

31-40    FB  -  Burette  calibration 

factor  ccs/division 

CHART  INPUT  (3F10) 

Columns   1-10   DIVST  -  Number  of  divisions  along  chart 

11-20    DIVLD  -  Number  of  load-cell  divisions 

21-30    DIVBU  -  Number  of  burette  divisions 

31-40    ICODE  =  0    test  factors  do  not  change 

=  1    after  this  point  new  test 

factors 

=  999  last  data  point 


There  will  be  a  type  iv.  card  for  each  data  point.   If 
ICODE  =  1  the  next  data  must  be  type  iii,  the  new  test 
factors . 
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C.3  Example  I' rob  1  em  for  Tri axial  Test  Program 
A  fictitious  triaxial  test,  Number  99,  has  been  run  and 
Table  C.l  lists  the  output  obtained  from  the  Sandborn 
Recorder . 

The  test  parameters  were 

Triaxial  machine  loading  rate   4.0  x  10" 5   ins/min. 

Sandborn  chart  speed  6.9  x  10"3  divs/min. 

Load-el]  calibration        0.8  and  2.0    lbs/div. 

Burette  calibration  0.1       ccs/div. 

and  it  was  desired  to  have  results  printed  out  in  both  tabu- 
lar form  and  on  all  eight  possible  plots.   Table  C.2  gives  a 
listing  of  the  sixteen  data  input  cards.   Note  that  during 
the  test  it  was  necessary  to  alter  the  load  output  scale  by 
switching  the  attenuator  to  a  higher  value.   It  was  therefore 
necessary  to  input  a  new  test  parameter  card,  the  ninth  card 
in  the  listing.   The  resulting  stresses  and  strains  arc 
listed  in  Table  C.3  in  tabular  form  while  Figures  C.2  and 
C.3  show  the  CALCOMP  plots  for  this  test. 
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TABLE  C.l 
Sandborn  Data  for  Example  Triaxial  Test  Reduction  Program 

Number  of  divisions   Number  of  load-     Number  of  volume 

along  chart  cell  divisions  burette  divisions 

0  0.0  0 

10  19.6  6 

20  30.1  10 

30  36.8  14 

40  41.8  17.5 

50  46.0  20.5 

CO  19.9  23 

70  21.3  25 

80  22.2  27 

90  23.2  28 

100  24.3  29 

110  25.4  29.6 

120  26.0  30 
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TAB LB  C.2 

Input  Data  Cards  for  Example  Triaxial  Test  Reduction  Program 

Card  No.  Card  Columns 

12 34 56 789 01 234 5678901234 S678901234 567890 

1  TRIAXIAL  TEST  NUMBER  99 

2  12345678 

3  0.0069  0.8  0.00004  0.1 

4  10.0  19.6  6.0  0 

5  20.0  30.1  10.0  0 

6  30.0  36.8  14.0  0 

7  40.0  41.8  17.5  0 

8  50.0  46.0  20.5  1 

9  0.0069  2.0  0.00004  0.1 


10  60.0  19.9  23.0  0 

11  70.0  21.3  25.0  0 

12  80.0  22.2  27.0  0 

13  90.0  23.2  28.0  0 

14  100.0  24.3  29.0  0 

15  110.0  25.4  29.6  0 

16  120.0  26.0  30.0  999 


0.8 

0.00004 

19.6 

6.0 

30.1 

10.0 

36.8 

14.0 

41.8 

17.5 

46.0 

20.5 

2.0 

0.00004 

19.9 

2  3.0 

21.3 

25.0 

22.2 

27.0 

23.2 

28.0 

24.3 

29.0 

25.4 

29.6 

26.0 

30.0 

I'M-, 


TABLE  C.3 
Output  from  Example  Triaxial  Test  Reduction  Program 

Triaxial  test  number  99 
Initial  cross  section  area  =  1.53918 
Initial  sample  volume      =  73.146 
Initial  sample  length      =  2.9  0 

Chart  speed  factor  (no.  divs/minute)    =  .007 

Load  factor  (no.  pounds/division)       =  .8 

Strain  rate  (no.  inches/minute)        =  .00004000 

Volumetric  strain  factor  (ccs/division) =  .10000 

Strain   Load   Buret   Time    Load   Disp.     Area    Del  Vol 
Divs    Divs   Divs    Hr     Lbs     Ins    Correct    ccs 

10.0  19.6  6.0  24.15  15.680  .05797  1.57058  .6000 

20.0  30.1  10.0  48.31  24.080  .11594  1.60328  1.0000 

30.0  36.8  14.0  72.46  29.440  .17391  1.63737  1.4000 

40.0  41.8  17.5  96.62  33.440  .23188  1.67295  1.7500 

50.0  46.0  20.5  120.77  36.800  .28986  1.71011  2.0500 

Chart  speed  factor  (no.  divs/minute)    =   .007 

Load  factor  (no.  pounds/division)      =  2.0 
Strain  rate  (no.  inches/minute)         =   .00004000 
Volumetric  strain  factor  (ccs/division) =   .10000 

60.0   19.9   23.0   144.93   39.800   .34783   1.74895   2.3000 
70.0   21.3   25.0   169.08   42.600   .40580   1.78960   2.5000 


80.0   22, 

.2   27.0 

193.24   44, 

.400   .46377 

1.83218 

2.7000 

90.0   23, 

.2   28.0 

217.39   46 

.400   .52174 

1.87684 

2.8000 

100.0   24. 

,3   29.0 

241.55   48 

.600   .57971 

1.92373 

2.9000 

110.0   25. 

,4   29.6 

265.70   50 

.800   .63768 

1.97303 

2.9600 

120.0   26, 

,0   30.0 

289.86   52 

.000   .69565 

2.02492 

3.0000 

Strain 

Stress 

Strain 

Strain  Octahedral 

Stresses 

Axial 

Axial 

Vol 

Dist 

Normal 

Shear 

.019990 

9.984 

.008203 

.017256 

3.328 

4.706 

.039980 

15.019 

.013671 

.035423 

5.006 

7.080 

.059970 

17.980 

.019140 

.053590 

5.993 

8.476 

.079960 

19.989 

.023925 

.071985 

6.663 

9.423 

.099950 

21.519 

.028026 

.090608 

7.173 

10.144 

.119940 

22.757 

.031444 

.109459 

7.586 

10.728 

.139930 

23.804 

.034178 

.128537 

7.935 

11.  221 

.159920 

24.233 

.036913 

.147616 

8.078 

11.424 

.179910 

24.722 

.038280 

.167150 

8.241 

11 .654 

.199900 

25.263 

.039647 

.186684 

8.421 

11 .909 

.219890 

25.747 

.040467 

.206401 

8.582 

12.137 

.239880 

25.680 

.041014 

.226209 

8.560 

12.106 
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L  I  ST  ING    OK    i'KuGhAm 


REAL    LOAD,    Lr.NUlH 

0  1  •'•^  «i>  I  j..    'j  I  vSI  (  5 'J  J)  .Dl  VlDI  Suili  .ul  v  duI  b  uj  )  ,  STRAIN!  500  )  • 

_1STjRcSS<  bJiO)  .DiLVJsOO)  ibHKS.IJt(.5O0_Lir.<>.00J  »u(  bOOJ 

DIM f  M 3 i  0*N    X Y  S  0 2 Y7 1 (50  2 I 
DIMFNSlON    TITLE! 3) 

call  plots _ 

RhADIbtl i  tTTLl 
FORMAT I 3A] 0 ) 
NPR1NT  -U 


CmOICF  OK  PLOTS 


1      IN    COLUMN     1 


»  IN    COLUMN 

t  IN_rOl.UMN 

■*  in  column 

>  In  column 

j  Hi  Cv«...  •:,. 

/  IN  COLUMN 

J  IN  COLUMN 


hlu I  Of 

r-LOl  Of 

>-LOT_  OF 

Pi  OT  OK 

PLOT  Or 

t-LuT  UK 

PLC  I  Or 


TAU  OCT  -  SIC  OCT 
jCI 


-  UhV  siKAIN 

-  vol  sikain.  _.. 

-  axi/-l  sixain 
li  -  axu.l  strain 

<jl    jIhAIn    -    j  111    OL  I     -    ^  1  ij    OC.I 
AXIAL     STKESS     -    AA|/»L     SlrfAlN 


TAU 

I  At  UC.I 

I  AU  ot I 

SlU 


PLOI  OK  VOL  olkAlN  -  Ut.V  SIKAIN 


Rc.A(.,<S  .2  I   1PLUT 
I OR^Al  I  J  8  I 


iPLf     UIMcNSlUNS 


__A".EA   .-_l»b391C 

VOL     *      /3.140 

LEN&Tii    =    ?.y 

WkIII.lt.  31     TilLL,     AREAi    VOL.     LtNC 

ro«"AT      I 1M1 » ///     7X.3A10.//// 

1  3?M0     INITIAL     LkGjS     StLtlON    ARtA 

2  32  MO     INITIAL     S«wPLl_VJLUMjl 

~J  ~32riO~"iNI  I  IAL   VmViKLc    LLNvJrii 


.K  l.b 
j>r-6._3  . 


_.  lLOOL=i 

1=0 
1  =  1  +  1 
[F  ( .IC.pI/E.N.Ei.n GO    TO    7 

<EA0     IN     TEST     PARAMEIKRS 


;L        -       Cn/.,<!      Si-lEU    KACluK  (     NO.     U1VS     /     MINUTE      I 

LO/'U     Fdflv/k  (     1,0.     POUNDS     /     DIVISION     ) 


.-V 


SIR..  IN  RATr._ 
VOLU'lt  IK  IL  ~S  i  HA "l  N 


iCn-.b_  /_MiNolL_ 
AlIGk  "'"(CCS  /  DIVISION  I 


Rt AOI5.5I   FCi  r L .  FT.  r  V 

FORMAT  ("<il  10.0  j 

WR I  IE  I  6 .0  I    KC.  FL,  M.  FV 

FORMAT  (  /' CHAKtT  SPtfcD  KACToK  j_ 

2/'       LOAU  FACTuK  (  no.  POUNDS  /  1 

3/'        STRAIN  KATE   I  NO.   INCHES  /  » 


UIVS  /  MINUTE  I 


1 V  I S 1  ON  I 

1NUTE  I 
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'•/'  vulu;  •  t.  1 1:  1l    s'r/.l.   ^A^.Ion    i    uj    /    DlVlaluN    )      =    '.ti.j    /) 

i'.a'jc.o/      i/ivornii   l/Ivluiiii   uivl-uiii.    icuiic 

fuvA  i      (  3M0.O.i.i.l_tJ 

L'),.(.-tM  VI  i>(  i  i«r  L 
I  I  Mi  -.i/l  V.SI  I  1  >/FC 

o;  •,!>=■  T  [v-  »k  t 

S  r-'A  INt   1  )  -!Ji  jf/LtWj  i.( 
A»FACN-AkIlA/(  1  .J-MKmIN<  I  )  I 

Sllvr<.',(   I  )  =[  (JA  J /ARE  AC.P 

I  I  Mill'/'    1  :  N'h  /h'J  •  J 
DFl  I  AV-DI  ViVJI   I   )  *fV 

ALNGTnf  Lt  NuTH    -    [,iSP        . 

VOl.rAl!:.  Ai.P.*ALNvjfM*i6.J6706<» 
ULLV { i I =Ut I    [ Atf/VCL 

RMHITN  I    I   )  -  S  I  I-  •,  i_t|  (   I   )_-.;lL7  I  I   I  /3  .u 

M     I    )    --D   I   r  |_V-.  <",    V/  i'.O 

SO-Ji.P.  I  I  ^.;j| 

01    i    |  =  l  'jU/3.0  )  •  ..  I  ^LiSI    II 

If     (Ni-KlNT.Nr  ,u|     00    10    lu 

W«|   !>■  ((,  ,V) 

I-',  rfl  ...I       I   ;n  1  .  ■,*  .  •_.',  I  rt  ;  ',  •  i*X_»  •  Lw>'u'_  •  :)  X  .  'pyHcJJ  .•  ->*■  I  '  T  I  r-'.c  '  ioX.  '  LuAu'i  _. 

i(./,|,jI..I'.'ii».i.'u.l/'.»(i'ocl    VuL  '  /  /X»  ■  01  Vb>  i^i  'Ul  Vi'  >hX  >  'wl  Vo*  «t>X 
/  >  ■  MM  'nXi'LUri'ioXi'lNs'ioKi'  LOKntC  r  •  ijX»'Cti'/) 

Nf'l?  |  NT-bC 

1,1'KKil  -*4I-K  INl'-'i 

Vi"  Ml    I  -  .  11  I  01  Vi.1  (  I  I  .  0  I  VLDl  !  I  .01  VOUI  i  I  »  T  IMluK.LOAD.DI  SP» 

1  /■  I-'  ,xr;<  ,  rj|  i    r  AV 

Hji-'ma't    "u  l?'.  J  .2?  if.  li'M'u.iVr  l»,"7i.2Fll.!>Tf:iO.<il 

I  I      I  I  CWS.it  .Nt  .999  I        GO     lo    <. 

t|.i-.nT-I 

rn:-  i  »j  r  =*ij 

O'J     )  '        ]  ~  1  iNWM.f  T 

1'      l.-|Ui<  IM  .Nt  .01 GO    IP__l* 

«  H  I  ! F 1 0  .  Ill 

liyr..','.l      I  1M1  ,8X.  •  STr<A  |;,<  ,  7X»  •  blKFbS'  t7X»  "SlriAIiN"  »8X  ,  •  bIKAlNI  ,6X< 
P  OCT  A»if  ORAL      '.THE  SLL  b'  /  rH  . '  Ax  |  Al  •  ,  (JX  .  '  Ax  !«!_•  i  VA.  '  «UL  '  i  lu*.  '  u  l  jI.'.i.BX 
2 .  '.'JOP'-'.-'-L'  •  iX«  ■  SHtAK"  /  j 

NPRlfiT-50 

_.fir'»l'-,J.=_'iPHjNj;rJ 

«,'*'  I  ri    (  6  .  1  'j  )  ST K A  INI   I  >  •  STkE4SI  I  I  .DtLVl  i  I  .SMKoll,  I  I  )  ,P(  I  )  ,0<  I  ) 

fOkvM     (FlA.6.F12.3.2f  lA.6t2F12.31 

CONTINUE , . :_..    

IF  (  H'LOT  .  tO.O  )  GO    TO    99 


M=NUMPT 

IF  I  II'LOl  .L  r. 2999999)  GO     TO     3/ 


PM.OT     OF      IAU    OCT 


f)0    1*.        1  =  1  .M 

xi  1 1  =p~i  i  r 

Y  (  ]  )  -  0  (  1  I 

xi m*-:  )=o.o 

X(M+2>=1.5 

Y I M  t ]  1=0.0 
Y(  M+2  1^_3.0 

xvn  im = xTk  i  1 1 

DX=X(M*2) 
YMIX)  =  Y(M+1  ) 
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DY  =  Y(M>-21 

CALI     AXIS     (0.0  ,0.0  ,?'.i>OCTAmLl>hAl    NonMAl    5TULi>S>,-2<i. /.  73  »U«UiXMIN, 

1L)X,UI      _     _     

CALL    ax  I  i     IO>u<U>Ui?3i10Cf  AHtl/HML    SHI-.AK    STRESS  1 23  »5.00iVO»0  •  YM1N»UY 
1  .01 

(.ALL    SYMIJGL  I  1.0>5..(J.O.  lOt.TI  I  LX.iu.Oi. 30) . . 

CALI     liNl     (  X.V..V,  1  ,-1,3) 
call    PLOT  (  lu.Ui'J, 0,-3) 
37_     IF  (JPL-.)T  .L  T.    J'fil'Ji)  CQ„IP_3ti 

PLOT     01       1AU     OCT     -     DlV     SIXAIN 

bo"  17  "i=  i  Tm   ~ 

XI  I  ) =5HRSTN<  I  ) 
_J7__Y<  IJ=<M  1  I 

x(M+i"f='6.o"' 

X(M+?)=0.03 

YIM-M )=0.0_ 

Y(M+2|  =  3".'0 
XM  J  n  =  X 1  M+ 1 ) 

YI1lf;  =  Y  (  M  +  1  ) 

DY  =  VI''>2I 

CALL    »x  I  :.  I  1.0.0. 'J .  1  7,,t  r  V  I /•,  I  C.<  K    S  I  KA  I  ft  ,  -17  ,  7  .  7b  ,  0  •  0  ,XMI  N.tOX  ,  0  I 

c-LL    <-x  <  ■>    (0.w«u.0t23i:Gi.TAnLuRi«L    bML<.r(    S  1  k:_  ss  •  ..  3  ,  5  .Co  >>0  .0«  YMI  N, 
1UY.-! I 

cali    s«"nUL  ( i .  o.«>.o.>  .  iu.  r ;  ilk  ,o..o«  io> 

CALL    L  llif     (  x.  r  ,,v,,"]  ,-1,3) 
C  ALL     PLO  I  I  lQ.O>u.0,-3) 

3h       If  I  I'  I  JT  .L  I  .       <.'V9991  CO    10    SV__ 

JL0T     Cf      'All    OC  I     -     VOL     SIHAIN 


DO    IB       I  *  1  ,M 
XI   I  > - Ll  LVII  ) 
1£_      Yl  |  )-'.!   II      

X  I  "l  *  1  1^0.0 

XI  t-U?  1  =  0.006 
_  Y(  'ft  1  )  -0,0 

Y(M+2)=3.0 
XMlN*Xf M+l ) 

DX=X(K+2)  ____    

VM I N  =  < I M  + 1 ) 
DY=Y<M+2) 

CALL     AX  1  S     10.0,0.0,  1  7MV0LUMETRIt_ST(<AIrt  ,_-  1  7j>  8.  CO  ,  0  .  0.  XM  l_N«_DX  .  OJ 

CALL    AX  I  b     <0.u,u.O,2AiiOCTRAHtDKAL    SrttAK    STRt  Si ,  Z<> ,  b  .  00  t  90 .  0  ,  YMl  N  , 
10V.-1 I 

CALL    SYMbOL(.i'0,5.0.0..1C,Tl  TlE.O.C.30) , . 

CALL  LINE  IX.Y.M, 1 ,-1 ,3| 
CALL  PLOT  I  10.0 ,0.0,-3  1 
39    If  I  |PL:jT  .L  T.     b999)     GO  TO  4U  


'LOT  OF  TAU  OCT  -  AXIAL  STRAIN 

DO  19   1=1 ,M  "~ 
X( I )>STRAIN( II 
Y(  I  |=0I  1  I 
— XTM+n  -=070  _ 

XIM+2)=0.'03 
YIM+1 )=0.0 
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Yl>"  +  ?  l^'i.O 
XM1N=XI M+] ) 

ox-y  (M »  pi 

VMlN  =  V ( M+  1  ) 
DY=Y(V  +  2  ) 
<  ALL     *<ll.l  ';.  I 
C  A  L I     AXIS     (  ( ) . 
1  ."I 
_    LAi.t     S/MuULI   | 
CALL     I    INt      I  X , 
'-ALL     PLOT  I  10  • 
If  I  |  PLOT.  I.  I  . 


(AXIAL     S  I  RmI  Mi-l*  .  /.  fS.O.U.XMl  lituttQ  ) 
SnuCTAnEuhAL    SHc  Ak    '.InLiii^lii.OUiVu.JiYMIhtUI 


)  t  5.0  i  U,  IOj.1  I  1  Lr. .  u..0«  JOJ 

.  M  .  1  .  -  1  ,  i  ) 
.  u  .  0  .  -  3  I 

6'JVI  GO    10    Al 


'LOT     Of      3IL>    Ol  1     -     AXl/.L     SIRrtiN 


DO    ?0       1-1 ,M 

X  I  I)  -:.  TRAI  Nl  I  ) 
Y(  | _|ip|  I  ) 

X  (  M  *  |    |   =  (j  .  0 

X(M-*2)  =  0.Ci3 
VfM*l ) :0.u 

Y(  »■>  /  I  -<?.U 
X  V  I N    XIM«II 

ux  =  x  i M.  •_  2_i 

Y  M  |  N  -  Y  (  M  f  1  ) 
UY  =  Y1.V*2  ) 


CALL     AXI5_!_0.Q  •0_._0_«  1  ?MAX  IAL_  5  I  K«  I  N  .  -  1  2  ,J  .  15  ,  0  .  0  .  XM  I  N  ..  fjX .  0 .).  _. 

CALL    AXIS     (O.'^t  v  ."  .^nuCl  AntL>.s«L    Nok.-.Al     Ji»til .  ft  t  5  .  00  »  VO.  0  t  YMlNi 
) OY.-l  I 

_       CALL    SrMUOL(  l.Ci.5.C.O.)0»TIJLe.0..p.30l 

CALI      L  I.V       (X,  Y.M,  1  ,-1,3) 
LA^L     k^Oi I |0,u . j.u,- j ) 

>,1        IFIIOLOT.tT.  .79) <>0._T0._A2. 


PLOT     OF     VOL     ',  TRAIN     -     S  I  0    OCT 


DO  ?  1   I  ■=  ]  ,  M 
XI  I  ) =P I  1  I 

_2JL.  Yl  |  )  .  Ur.  LV(J) 

X1M+1 1=0.0 
X(W*?)=2.0 

Y(M»1)=0.0  _ 

""  Y(M+2)-0.01 

XM1N=X(M*1 ) 

DX  =  x  ('■■*?  I 

""  "yMIN='yIMM  ) 

DY=Y(M+2 I 

CALL  AX  I  ',  (  C.O.  0.0,  2AHOCIAMCDHAL  NOKSA 
1,01 
Call  Axis  (0.0,0.0, 17HV0LUMETRU  ST  H,\  I  r< ,  1  7  ,  5  .  00  .  90  .0  •  YM  I  N  ,UY  ,  -1  ) 
CALL  S  'r-.'.uL  (  1  .O.S. 0.0. 10.1  1  TLE  .0.0. i0| 

Call  l Inf. "  I  x  .  v ',m", "i","-i  .  3  ) 

CALL  PLOT ( 10.0,0.0.-3) 
Ai    IK11PL01.LT.        9)     GO  TO  9t) 


,TkESS»-2<..  /.7S.0.0.XM1M.DX 


5L0T  OF  AXIAL  STkcSS  -  AXIAL  STRAIN 


DO  ??    1  =  1  ,M 
XI  1  )=S TRAIN!  I  I 
2  c       YlI)=STRtSS(l) 
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X(  •*+ 

XI  M* 

v  i  u* 

Y  [  V. 

XVI  N 
DX  =  X 

r  k  l  n 

OY=Y 

CALL 
"CALL 
CALL 
CALL 
"CALL 


1  t=U. 
?l=0. 

1  >  =  n. 

Zf=&. 
=  X  (  M  - 
<M+2) 
=  Y  IM» 
<>•'.♦?> 
AXIS 
AxT'S 
SYRs 


0.0.0.0 t ]  ?HAX1  AL  £.1'*A1N  t_-l*l  »  '. 

c.i.o . j. "■.'/."] ?maa j al  MKtisrr<r>i. 

LI  l.l.5.O.0.10iTlTLci0.0t3O) 
LiKii     I  X.  Y. /. ]  ,-1  ,3) 
PLOT  (  IO'.O'iO.O  ,-3"| 


i  <  X  M  I  N  .  D  X  .  0  )  _ 

rdiYMiNior  i-i  r 


->LOT     OF     VOL     STRAW,     -     DEV    STRAIN 


DO    ?3       l=l.« 
XI 1 |=SmHSTN( 1) 
Y(  I  )  =lJtC  I  VI  I  )" 
xi  "i-i  )  =0.0 
X(M+7)=o.03 

V(  M»l   l-U.IJ 

Y( »*s ) =0.0] 
XMIN»X  I  MM  | 
!ii  =  «l">2l 
YMlN'Yl M > 1  I 
0Y»VIM+?I 

Call  Vix  I  j(  i.r.j. 

CALL  AXI'.  10. vj 

CALL  SYMBOL  I  1  • 

(.ALL  Li  m   I  >   .  t 

CALL  f'L'jT  I  lC.vJ 

CALL  l-l-O  I  (Ot'J. 
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APPENDIX  D 
THE  FINITE  ELEMENT  CONSOLIDATION  PROGRAM 
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D.   THE  FINITE  ELEMENT  CONSOLIDATION  PROGRAM 

D . 1   Input  Format  for  Consolidation  Program 

i.   IDENTIFICATION  CARD   (18A4) 
Columns   1-72   HED  (18) 

ii.   CONTROL  CARD   (815,  F10.0,  215,  F10.0) 

Columns   1-5   NUMNP  -  Number  of  nodal  points 
6-10   NUMEL  -  Number  of  elements 
11-15   NUMAT  -  Number  of  different  materials 
16-20   NUMPC  -  Number  of  boundary  loading  cards 
21-25   NUMSP  -  Number  of  elements  with  speci- 
fied pore  pressures 
26-30   NUMBC  -  Number  of  flow  boundary  cards 
31-35   NDECR  -  Code  for  degree  of  consolida- 
tion calculations 
36-40  NFPERC  -  Percentage  desired 
41-50   FTIME  -  Final  consolidation  time 
51-55  NLOADC  -  Code  for  time-dependent  loading 
56-60   NPLT  -  Code  for  plotting 
61-70   DTIME  -  Initial  consolidation  time 
increment 
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iii.   MATERIAL  PROPERTY  CARDS   (15,  5X,  4F10.0) 
Columns   1-5   MTYPE  -  Material  number 
6-10  Blank 

11-20      E  -  Modulus  of  elasticity 
21-30    GNU  -  Poisson's  ratio 
31-40  COND,l  -  Horizontal  permeability/unit 

weight  of  water 
41-50  C0ND,2  -  Vertical  permeability/unit 
weight  of  water 
There  will  be  a  type  iii  card  for  each  material  type. 

iv.   NODAL  POINT  CARDS   (215,  4 PI 0.0) 

Column    1-5     N  -  Nodal  point  number 

6-10  CODE  -  Code  for  specified  forces  or 
displacements 
11-20    X  -  X-coordinate 
21-30    Y  -  Y-coordinate 

31-40   UX  -  Force  or  displacement  in  X-direction 
41-50    UY  -  Force  or  displacement  in  Y-direction 
There  will  be  one  card  for  each  nodal  point 

v.   ELEMENT  CARDS   (615) 

Column    1-5     M  -  Element  number 

6-10  IX, 1  -  Nodal  point  I 

11-15  IX, 2  -  Nodal  point  J 

16-20  IX, 3  -  Nodal  point  K 

21-25  IX, 4  -  Nodal  point  L 

26-30  IX, 5  -  Material  type 
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vi.   BOUNDARY  LOAD  CARDS   (215,  F10.0) 
Columns   1-5   IBC  -  Nodal  point  I 
6-10  JBC  -  Nodal  point  J 
11-20   PR  -  Normal  pressure 
There  will  be  NUMPC  cards  of  type  vi. 

vii.   ELEMENT  PORE  PRESSURE  CARDS   (215,  F10.0) 
Columns   1-5      L  -  Element  number 

6-10    NL  -  +ve  -  pressures  initially  only 

-ve  -  pressures  apply  throughout 
consolidation 
11-20  PRESS  -  Specified  pressure 
There  will  be  NUMSP  cards  of  type  vii 

viii.   FLOW  BOUNDARY  CARDS  (415,  F10.0) 

Columns   1-5     I  -  First  boundary  node 

6-10    J  -  Second  boundary  node 

11-15    K  -  Element  number  inside  boundary 

16-20  KODE  -  Code  -  flow  or  no  flow 

21-30   BPR  -  Pressure  outside  boundary 

There  will  be  NUMBC  cards  of  type  viii. 

ix.   DEGREE  CARDS 

x.   TIME  DEPENDENT  LOADING  CARDS 

xi.   PLOT  CARDS 
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